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INTRODUCTION 

The  longstanding  technological  demand  for  improved 
performance  in  materials  continues  unabated .  In  fact,  the 
promise  of  more  than  a  decade  in  filamentary  composite  ma¬ 
terials  has  generated  intense  interest.  Such  interest  is 
well  founded,  having  been  generated  by  the  possibility  of 
using  brittle  materials  with  high  modulus,  high  strength, 
but  low  density  in  composites  which  fail  in  a  non-cata- 
strophic  manner.  Such  fiber  reinforced  composite  materials 
offer  substantially  improved  performance  and  potentially 
lower  costs  for  aerospace  hardware. 

Fulfilling  the  promise  of  composite  materials,  however, 
requires  a  strong  technology  base.  NASA  and  AFOSR  have  re¬ 
alized  that  to  fully  exploit  composites  in  sophisticated 
aerospace  structures  the  technology  base  must  be  improved. 
This,  in  turn,  calls  for  expanding  fundamental  knowledge 
and  the  means  by  which  it  can  be  successfully  applied  in 
design  and  manufacture.  It  also  calls  for  expanding  the 
body  of  engineers  and  scientists  competent  in  these  areas. 

As  part  of  their  approach  to  accomplishing  this,  NASA  and 
AFOSR  have  funded  the  current  composites  program  at  Rens¬ 
selaer.  The  purpose  of  the  RPI  composites  program  is  to 
develop  advanced  technology  in  the  areas  of  physical  prop¬ 
erties,  structural  concepts  and  analysis,  manufacturing, 
reliability  and  life  prediction.  Concommitant  goals  are 

Preceding  page  Wank 
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to  educate  engineers  to  design  and  use  composite  materials 
as  normal  or  conventional  materials.  A  multifaceted  pro¬ 
gram  has  been  instituted  to  achieve  these  objectives.  The 
major  elements  of  the  program  are: 

1.  CAPCO^g  (Composite  Aircraft  P^rogram  Component)  . 
CAPCOMP  is  primarily  a  graduate  level  project  being  con¬ 
ducted  in  parallel  with  a  composite  structures  program 
sponsored  by  NASA  and  performed  by  a  private,  aerospace 
manufacturing  contractor.  The  first  component  redesign  is 
being  done  in  conjunction  with  the  Boeing  Commercial  Air¬ 
plane  Company.  The  main  spar/rib  region  on  the  Boeing  727 
elevator,  near  its  actuator  attachment  point,  has  been  se¬ 
lected,  with  Boeing's  advice  and  the  concurrence  of  NASA/ 
AFOSR,  for  study  in  CAPCOtiP .  The  magnitude  of  the  project 
-  studying,  designing,  fabricating  and  testing  a  relatively 
small  but  highly  stressed  region  on  the  elevator  -  is  both 
consistent  v.’ith  Rensselaer's  capabilities  and  a  significant 
challenge.  The  selection  of  a  portion  of  a  full-scale 
flight  hardware  structure  assures  relevance  to  this  pro¬ 
ject's  direction. 

Visits  to  Boeing  were  conducted  in  the  Fall  of  1978 
by  Professor  Hoff  and  several  of  his  students,  and  the  first 
serious  design  work  began  shortly  thereafter.  Two  alterna¬ 
tive  designs  were  pursued  to  the  point  of  preliminary 
analysis  and  testing;  further  progress  is  reported  in  Part 
I. 
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2.  CAPGLIDE  (Composite  Aircraft  Program  Glider) . 

This  undergraduate  demonstration  project  is  to  design,  fab¬ 
ricate  and  test  a  foot-launched,  ultralight  glider  using 
composite  structures.  A  flight  vehicle  was  selected  to 
maximize  student  interest  and  to  provide  the  students  with 
a  broad-based  engineering  experience.  For  those  students 
continuing  with  graduate  work  at  PJ>I,  CAPGLIDE  is  intended 
to  provide  natural  progression  to  CAPCOMP.  The  progress 
on  the  CAPGLIDE  project  to  date  has  been  very  good.  Seven 
professors  and  approximately  30  students  were  actively  en¬ 
gaged  in  the  project  during  the  beginning  of  this  period; 
that  is,  at  the  end  of  the  spring  semester.  A  description 
of  the  status  of  the  work  performed  under  CAPGLIDE  to  the 
end  of  this  reporting  period  is  given  in  Part  II. 

3.  COMP AD  (Computer  Aided  Design) .  A  major  thrust 

of  the  composites  program  is  to  develop  effective  and  effi¬ 
cient  tools  for  the  analysis  and  design  of  composite  struc¬ 
tures.  Rensselaer  and  NASA  Langley  have  jointly  implemented 
the  use  of  the  SPAR  code  on  minicomputers,  and  the  work  at 
Rensselaer  has  made  "virtual  memory"  available  to  those 
using  SPAR.  Recently  attention  has  turned  to  preprocessor 
developments;  details  are  reported  in  Part  III. 

4.  Composites  Fabrication  and  Test  Facility.  Struc¬ 
tural  design  engineers,  educated  only  by  course  work  and 
design  projects  limited  to  paper,  often  fail  to  sense  or 
appreciate  problems  involved  in  fabrication.  The  actual 
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fabrication  and  testing  of  composite  structural  components 
provides  this  training  and  the  final  validation  for  the  de~ 
signs  in  our  CAP  projects.  PJ>I's  Composites  Fabrication 
and  Test  Facility  is  located  in  the  laboratory  and  high  bay 
areas  of  the  Jonsson  Engineering  Center.  Equipment  is 
available  for  compression  molding  parts  as  large  as  19"  x 
19"  and  vacuum  bagging  parts  of  much  larger  size.  Panels 
approximately  4'  x  20'  have  been  made  by  vacuum  bagging. 
Specifics  are  given  in  Part  II,  CAPGLIDE.  A  pressure  vessel 
for  small  parts  and  spars  has  been  designed  and  built.  A 
second  one  capable  of  higher  pressures  is  also  planned;  more 
complete  details  are  reported  in  Part  I,  CAPCOMP . 

Our  capabilities  and  efficiency  in  fabrication  and 
■^^03'ting  were  increased  substantially  during  the  reporting 
period  by  the  addition  of  a  Manager  and  Master  Technician 
for  the  Composites  Laboratory.  This  position  is  filled  by 
Volker  Paedelt  whose  background  includes  almost  twenty  years 
at  the  Vereinigte  Deutscne  Metallwerke,  Frankfurt,  Germany; 
Simmons  Fastener,  Menands,  N.  Y.  and  the  Physics  Department 
at  Rensselaer  Polytechnic  Institute 

5.  INSURE  ( Innovative  and  Supporting  ^search)  .  The 
criteria  for  selection  of  research  projects  to  be  conducted 
under  this  program  are  (a)  that  they  must  anticipate  criti¬ 
cal  problem  areas  which  may  occur  in  the  CAP  or  NASA/AFOSR 
programs  or  (b)  that  solutions  to  existing  problems  are  not 
yet  satisfactorily  in  hand.  During  the  reporting  period 
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five  programs  were  funded,  Res-ults  from- the  ongoing  projects 
are  reported  in  Part  IV. 

6,  Curriculum  Revisions,  The  goal  of  educating  engi“ 
neers  to  think  of  composites  as  normal  or  conventional  mate¬ 
rials  has  required  changes  in  curriculum.  Since  the  initi¬ 
ation  of  this  program,  almost  all  Rensselaer  engineers  take 
introductory  courses  which  incorporate  the  concepts  of  ani¬ 
sotropy  and  composite  materials.  In  addition,  six  special¬ 
ized  courses  in  composites  have  been  offered  during  the  past 
three  years  to  develop  those  special  skills  required  of 
students  involved  in  the  composites  program.  A  new  course 
v;as  introduced  in  the  Fall  of  '78  semester  on  composite 
design  and  analysis  using  programmable  hand  calculators,  a 
central  mini  and  full  frame  computers.  A  new  graduate  level 
advanced  topics  course  with  the  title  "Advanced  Finite  Ele¬ 
ments"  was  offered  for  the  first  time  in  September  1979. 

The  additions  of  the  SAP  and  SPAR  computer  codes  and 
the  growing  availability  of  interactive  computer  graphics 
under  our  COMPAD  program  element  have  reached  the  point 
where  our  engineering  students  are  using  these  facilities 
as  everyday  working  tools  for  design,  analysis  and  visuali¬ 
zation  purposes.  We  have  thus  achieved  one  of  the  principal 
goals  of  the  curriculum  development  activities, 

7 .  Technical  Interchange. 

a)  Dr.  Kiyoshi  Kenmochi  returned  to  the  Industrial  Products 
Research  Institute  of  Japan's  Industrial  Science  &  Technology 
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Agency  in  August.  He  was  most  'helpful  during  his  appoint¬ 
ment  as  Research  Associate  under  the  program,  having  set  up 
our  photoelastic  stress  analysis  facilities  and  assisting 
us  in  conducting  the  mechanical  attachment  research  reported 
in  our  last  report  (July,  1979)  and  herein, 
b)  Technical  Meetings:  Technical  meetings,  on-  and  off- 
campus,  provide  important  opportunities  for  interchange  of 
technical  information.  Because  of  the  large  number  of  com¬ 
posites  meetings,  a  central  catalog  with  all  upcoming  meet¬ 
ings  is  being  maintained  and  distributed  periodically.  In 
this  way  we  help  assure  that  a  Rensselaer  staff  member  will 
participate  in  important  meetings.  The  calendar  for  this 
reporting  period  is  shown  in  Table  I.  Meetings  attended  by 
RPI  composites  program  faculty/staff  during  the  reporting 
period  are  shown  in  Table  II.  Some  meetings  particularly 
relevant  to  composites,  held  on-campus  with  off-campus 
speakers,  are  listed  in  Table  III.  A  list  of  composite-re¬ 
lated  visits  to  relevant  organizations  by  RPI  faculty/staf f / 
students,  with  the  purpose  of  each  visit  outlined,  is  pre¬ 
sented  in  Table  IV.  Regular  readers  of  these  progress  re¬ 
ports  will  note  the  higher  level  of  such  activity  in  the 
months  from  September  to  April  -  partially  because  full  aca¬ 
demic  year  programs  are  in  session  -  than  in  those  reports 
covering  periods  including  the  summer. 

In  summary,  the  NASA/AFOSR  Composites  Aircraft  Program 
is  a  multi-faceted  program  whereby  aeronautical,  mechanical 
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TABLE  I 

CALENDAR  OF  COI4POSITE'S-RELATED  ' MEETINGS 
for  the  period  April  '79  to  September  '79 


1979 

4/4-6 

4/16-27 

5/7-10 

5/9-11 

5/14-17 

5/21-25 

5/21-25 

5/22-23 

5/23-25 

6/14-16 

6/18-20 

6/25-29 

6/25-29 

7/11-13 

8/6-10 

8/20-22 

9/3-7 

9/10-12 

9/10-13 

9/11-14 


20th  Structures,  Structural  Dynamics  &  Materials 
Conf.,  St.  Louis,  Mo.  Sponsored  by  AIAA/ASME/ASCE . 

Very  Large  Vehicle  Conference,  Arlington,  Va. 

Sponsored  by  AIAA. 

Design  Engineering  Conference  &  Show,  Chicago,  Ill. 
Sponsored  by  ASME. 

Materials  Conference,  Washington,  D.  C.  Sponsored 
by  ASME. 

Conference  on  Space  Manufacturing  Facilities,  Prince¬ 
ton,  N.  J.  Sponsored  by  Princeton  and  AIAA. 

Fracture  of  Composites,  San  Francisco,  Ca.  Sponsored 
by  Society  for  Experimental  Stress  Analysis. 

High  Modulus  Fibers  &  Their  Composites,  San  Fran¬ 
cisco,  Ca.  Sponsored  by  ASTM. 

Symposium  on  the  Fatigue  of  Composite  Materials, 

San  Francisco,  Ca.  Sponsored  by  ASTM. 

Advanced  Composites:  Design  and  Applications, 
Gaithersburg,  Md.  Sponsored  by  National  Bureau  of 
Standards,  NASA  and  MADC. 

Reinforced  Plastics/Composites,  Chicago,  Ill.  Spon¬ 
sored  by  Plastics  Seminars. 

Applied  Mechanics,  Fluids  Eng.  and  Bioeng.  Conference, 
Niagara  Falls,  N,  Y.  Sponsored  by  ASME  and  CSME. 

Pressure  Vessels  &  Piping  Conference,  San  Francisco, 
Ca.  Sponsored  by  ASME. 

Conference  on  Carbon,  University  Park,  Pa. 

Lighter -than-Air  Systems  Technology  Conference,  Palo 
Alto,  Ca.  Sponsored  by  AIAA. 

Sixth  Annual  Conference  on  Computer  Graphics  &  Inter¬ 
active  Techniques,  Chicago,  Ill. 

Aircraft  Systems  &  Technical  Meeting,  New  York,  N.  Y. 
Sponsored  by  AIAA. 

Practical  Applications  and  Design  of  Fibre  Composites, 
University  of  Cambridge,  England. 

Design  Engineering  Technical  Conference,  St.  Louis, 

Mo.  Sponsored  by  ASME. 

International  Conference  on  Acoustic  Emission,  Ana¬ 
heim,  Ca.  Sponsored  by  ASNT . 

International  Forum:  Composites,  Assemblage  &  Adhe¬ 
sion,  Lyon,  France. 
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TABLE  II 

COMPOS ITES-PvELATED  TECHNICAL  MEETINGS  ATTENDED  OFF-CAMPUS 
for  the  period  April  '79  to  September  '79 

1979 

5/21-26  ASTM  Symposium  on  Composite  Materials,  Visit  EPRI 
and  Lawrence  Livermore  Lab.  (Prof.  Stoloff ) . 

5/21-22  American  Helicopter  Society  Forum  (Prof.  Loewy) , 
Washington,  D.  C. 

5/21-24  ASTM  Fatigue  Conference  on  Composites  (Prof.  Krempl) , 
San  Francisco,  California. 

5/7-8  Conference  at  National  Bureau  of  Standards  on  Tech¬ 
nical  Cooperative  Center  for  Composites  (Prof.  Le- 
Maistre) ,  Washington,  D.  C. 

7/11-13  Acoustical  Society  of  America  Meeting,  (Prof.  Das 
and  C.  Lanzl) ,  Boston,  Massachusetts. 


TABLE  III 

COMPOSITES -RELATED  MEETINGS /TALKS  HELD  AT  RPI 


(April  '79  - 

September 

'79) 

Topic 

Date 

Speaker (s) 

Finite  Element  Analysis  of 
the  Post  Buckling  Behavior 
of  Curved,  Laminated,  Com¬ 
posite  Panels 

4/9/79 

N.  J.  Kudva,  Virgin¬ 
ia  Polytechnic  Inst. 

Structural  Design  for 
Longivity:  A  Professor's 

View 

4/24/79 

James  W.  Mar,  Massa¬ 
chusetts  Institute 
of  Technology 

Stress  Analysis  Around 
Macroscopic  Cracks  Under 
Creep  Conditions 

4/27/79 

Hermann  Reidel,  Max 
Planck  Institute  for 
Iron  Study 

The  Effect  of  Statistical 
Correlation  on  Weakest 

Link  Systems 

5/3/79 

R.  C.  Garson,  Univer¬ 
sity  of  Pittsburgh 

Aromatic  Polyamides:  De¬ 
velopment  and  Uses 

10/10/79 

Paul  W.  Morgan,  for¬ 
merly  Research  Asso¬ 
ciate,  E.  I.  duPont 
de  Nemours  and  Co. 

New  Methods  in  High  Tem¬ 
perature  Structural  De¬ 
sign 

10/11/79 

Donald  Griffin,  Wes- 
tinghouse  Corporation 

World  Trends  in  Computer- 
Aided  Manufacturing 

10/25/79 

M.  Eugene  Merchant, 
Director  of  Research 
Planning,  Cincinnati 
Milacron,  Inc. 
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TABLE  IV 


COMPOSITE-RELATED  VISITS  TO  RELEVANT  ORGANIZATIONS 


by 

RPI  Faculty/Staff/Students 

visited 

Date 

by  Prof (s) 

Purpose 

AFOSR: 

Dr.  L.  Kravitz, 
Director 

5/10/79 

R.  Loewy 

To  discuss  NDE, 
Fatigue  and 
Failure  Mechan¬ 
isms 

W .  P  .  A .  F .  B  .  ; 

Mr.  G.  Peterson 
and  personnel 
of  AFML  and 
AFFDL 

7/30/79 

J.  Diefendorf 

Lo  Feeser 

G.  Helwig 

R .  Loewy 

S.  Sternstein 

To  discuss  the 
NASA/AFOSR  Pro¬ 
gram  and  USAFML 
amd  USAFFDL  pro 
grams  in  compos 
ites 
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and  materials  engineers  must  interact  to"  achieve  its  goals. 
"Hard-nosed"  engineering  of  composite  aircraft  structures 
is  balanced  against  research  aimed  at  solving  present  and 
future  problems.  In  the  following  sections,  detailed  de¬ 
scriptions  of  the  CAPCOMP,  CAPGLIDE,  COMPAD  and  INSURE  pro¬ 
grams  are  presented. 


PART  I 


CAPCOMP  (Composite  Aircraft  Program  Component) 
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CAPCOMP  (Composite  Aircraft  Program  Component) 

(N.  Hoff,  R.  Loev;y) 

CAPCOMP  is  a  program  to  design  flight  critical  struc¬ 
tures  to  take  the  maximum  advantage  of  composite  materials. 

By  combining  the  efforts  of  experienced  faculty  with  bright 
and  well  trained  but  inexperienced  graduate  students  in  an 
environment  relatively  free  of  traditional  design  and  manu¬ 
facturing  processes,  we  hope  to  devise  new  and  hopefully 
useful  design  concepts. 

There  is  sufficient  information  available  today  to 
prove  that  many  structural  elements  can  be  built  lighter  of 
advanced  composites  than  of  metals.  But  if  such  elements 
have  to  be  joined  by  any  other  method  than  adhesive  bonding, 
difficulties  and  uncertainties  arise  which  can  be  eliminated 
only  through  conservative  designs  with  their  attendant  pen¬ 
alties  in  weight  or  by  extensive  and  expensive  programs  of 
"cut  and  try" .  This  stands  as  one  important  impediment  to 
full  adoption  of  composites  by  the  aerospace  industry. 

On  the  basis  of  these  considerations  Rensselaer  Poly¬ 
technic  Institute  began,  as  the  first  task  aimed  at  new 
structural  concepts,  the  design  using  composites  of  a  joint 
used  in  an  airplane  elevator.  To  make  the  design  realistic, 
an  existing  metal  airframe  component  was  chosen  for  redesign 
in  composites.  The  existing  design  chosen  was  that  of  the 
Boeing  727  elevator  actuator  attachment.  We  conceive  of 
this  work  as  carrying  forward  a  Structures  Demonstration 
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Program  using  the  joint  of  the  727  elevat:or,  paralleling 
that  of  NASA  and  its  aerospace  engineering  contractor,  the 
Boeing  Commercial  Airplane  Company.  Our  design,  fabrication 
and' test  effort  will  emphasize  design  ideas  specifically 
suited  to  advanced  composite  construction  for  the  purpose 
of  minimizing  the  v;eight  of  the  structure,  but  on  a  scale 
consistent  with  the  university  context  and  funding  level. 

The  staff  of  RPI  is  very  grateful  to  the  Boeing  Company  and 
its  engineers  for  their  wholehearted  support  of  the  work 
at  RPI. 

During  the  reporting  period,  the  two  different  designs 
reported  previously  (July,  1979)  as  suitable  for  replacing 
the  largely  metal  attachiaent  produced  by  Boeing  were  pur¬ 
sued.  The  first,  using  quasi-isotropic  graphite-epoxy 
laminates  (the  R.  W.  Berg  design) ,  was  chosen  for  initial 
fabrication  and  second-round  stress  analysis.  The  second 
(C.  Muser)  design,  in  which  a  deliberate  attempt  was  made 
to  use  uniaxial  graphite-epoxy  tape  to  as  great  an  extent 
as  possible,  was  suspended  in  favor  of  new  design  analyses 
and  test  planning  v/ith  the  objective  of  optimizing  the 
strength  of  the  movable,  heavily  loaded  bolted  joint  attach¬ 
ment  point. 

The  research  aimed  at  maximizing  the  load  carrying 
ability  of  the  attachment  of  the  rib  flanges  to  the  elevator 
skin  through  mechanical  (i.e.,  pin-type)  fasteners  also  con¬ 
tinued.  The  literature  and  theory  of  empty  and  loaded 
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circular  holes  in  composite  pla'tes  was  reviewed  as  part  of 

<9 

the  process  of  determining  the  most  efficient  arrangement 
of  fibers  around  the  hole, 

'  These  efforts  have  been  jointly  directed  by  Dr.  Nicholas 
J.  Hoff,  in  part-time  arrangement,  and  by  Dr.  Robert  G. 

Loewy.  They  are  described  further  in  the  follov.’ing  para¬ 
graphs  . 

1.  The  Elevator  and  Its  Attachment 

The  conventional  aluminum  alloy  elevator  of  the  Boeing 
727  is  shovm  in  the  upper  half  of  Figure  1.  The  lower  half 
of  the  figure  is  the  new  version  of  the  elevator  redesigned 
by  Boeing  in  graphite_epoxy ;  it  is  evident  from  the  pictures 
that  the  latter  is  composed  of  fewer  parts  than  the  former. 
However,  the  actuator  fitting  of  the  new  design  is  still 
manufactured  of  aluminum  alloy.  This  fitting  is  shown  in 
Figure  2.  The  fitting  is  attached  to  outboard  and  inboard 
portions  of  a  new  graphite-epoxy  spar  and  to  a  graphite- 
epoxy  nomex-honeycomb  rib  as  indicated  in  Figure  3, 

The  a tta^'hmeTrE'Ws' designed  by  Boeing  to  carry  loads 
up  to  19,000  lbs.  The  direction  of  the  load  varies  as  the 
elevator  rotates  over  an  angle  of  28  degrees  from  the  full- 
down  to  the  full-up  position. 


2.  Berg's  Design 


Berg's  design  is  shown  in  Figures  4  and  5.  The  first 


ALUMINUM  ACTUATOR  FITTING  (Boeing  Design) 


Figure 


COMPOSITE  ELEVATOR  ACTUATOR  ASSEMBLY  (Boeing  Design) 

Figure  3 
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ACTUATOR  LUG  TO  RIB  TRANSITION  -  TOP  VIEW 


Figure  4 


PROPOSED  SPAR  AND  RIB  ASSEMBLY  (Berg  Design) 
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of  these  figures  compares  the  Boeing  composite  design  with 
Berg's  design;  both  are  shown.  Figure  5  shows  how  the  spar 
and  new  actuator  rib  would  be  assembled  and  the  build-up  of 
thickness  in  the  lug  areas  anticipated,  as  needed,  to  provide 
the  required  bearing  strength.  Note  that  the  edges  of  some 
of  the  layers  are  bent  90  degrees  to  form  flanges  to  which 
the  upper  and  lower  coverplates  of  the  elevator  can  be  attach¬ 
ed.  Attachment  would  be  by  means  of  titanium  Hi-Loc  fasteners. 
The  right-hand  and  left-hand  graphite-epoxy  webs  are  stabiliz¬ 
ed  by  a  layer  of  nomex-honeycomb  between  them. 

The  finite-element  analysis  of  the  actuator  rib  v/eb 
(including  the  areas  of  the  lugs  with  the  bolt  holes)  re¬ 
ported  for  the  last  period  has  been  refined.  In  the  first 
computation,  the  wall  thickness  of  the  actuator  rib  web  was 
taken  as  constant  to  examine  the  stress  distribution  in  the 
area  of  the  holes  and  the  way  these  stresses  transfer  loads 
to  the  flanges  of  the  actuator  rib.  The  mathematical  model 
contained  482  nodes  and  a  total  of  800  quadrilateral  and 
triangular  elements  (see  Figure  6) ,  introduced  forces  simu¬ 
lating  bolt  loads  from  the  actuator  and  reaction  links  (see 
Figure  2)  and  a  distributed  shear  flow  along  the  upper  and 
lower  edges  of  the  actuator  rib  to  represent  the  reacting 
torque  resulting  from  aerodynamic  forces  on  the  elevator. 

In  the  current  reporting  period  the  analysis  was  re¬ 
peated  for  a  variable  thickness  web,  using  the  same  mathe¬ 
matical  models.  Several  iterations  were  made  examining 
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stress  levels  and  varying  thickness  distributions  until  a 
relatively  uniform  stress  level  design  resulted.  The  vari¬ 
able  thickness,  finite  element  model  is  shown  in  Figure  7. 
Graphic  comparisons  of  the  thickness  distribution  shown  in 
the  last  report  with  that  resulting  from  the  analysis  de¬ 
scribed  here  can  be  made  by  examining  Figures  8a  and  8b. 
Similarly,  comparison  of  the  stress  distribution  as  last 
reported  for  uniform  thickness  with  that  calculated  as  de¬ 
scribed  above  can  be  made  by  examining  Figures. 9a  and  9b. 

These  finite  element  analyses  results  revealed  two 
possible  unconservative  features  of  Berg's  initial  design, 
which  may  otherwise  have  gone  undetected.  First,  the  plate 
thickness  around  the  attachment  bolt  holes,  based  on  the 
Advanced  Composite  Design  Guide,  v/as  not  large  enough  to 
insure  the  desired  stress  levels.  The  more  detailed  model¬ 
ing  described  here  was  necessary  to  accurately  assess  these 
stresses . 

Second,  there  is  a  large  stress  concentration  at  the 
corner  just  aft  of  the  upper  bolt  hole.  The  Berg  thickness 
distribution  shown  in  Figure  7  had  not  yet  taken  account  of 
that,  even  after  the  first  finite  element  analysis  redesign. 

In  preparation  for  fabricating  the  first  full  elevator 
actuator  attachment  rib  to  be  undertaken,  step-by-step  plan¬ 
ning  of  the  fabrication  method  was  begun.  This  planning  is 
to  include  the  tooling,  the  (quasi-isotropic)  composite  ma¬ 
terial  lay-up  patterns  and  stacking  sequence,  the  curing 
cycles  and  estimates  of  the  kind  of  distortion  and  residual 
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VARIABLE  THICKNESS  FINITE  ELEMENT  MESH 


ACTUATOR  RIB  MODIFIED 
FOLLOWING  FINITE  ELEMENT 
ANALYSIS  (Berg  Design) 


MODIFIED  ACTUATOR  RIB  DESIGN 
INCORPORATING 
VARIABLE  THICKNESS  FINITE 
ELEMENT  ANALYSIS  RESULTS 
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Principle  Stress  Contours  for  Uniform  Thickness  of  0.375  Inches 
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stresses  that  would  be  expected  to  result  using  the  proposed 
procedure. 

In  the  first  stages  of  this  planning,  it  became  clear 
that  ease  of  operation  and  other  efficiencies  argue  strongly 
for  obtaining  a  small  autoclave,  just  large  enough  to  make 
the  kinds  of  parts  that  we  are  considering  under  CAPCOMP; 
i.e.,  between  one  to  two  feet  in  diameter  by  about  five  feet 
long  internally  and  capable  of  withstanding  pressures  in  the 
75  to  100  psi  range.  In  the  time  that  it  takes  to  bring  this 
new  facility  on  line,  experiments  with  the  proposed  flange- 
web  lay-ups  will  be  undertaken  using  vacuum  bagging  tech¬ 
niques,  to  investigate  distortion  problems. 

3.  Optimizing  Fiber  Orientations  in  The 
Vicinity  of  Heavily  Loaded  Joints 

It  has  often  been  noted  that  a  basic  advantage  of  fiber- 
reinforced  materials  is  the  possibility  of  "tailoring"  the 
material  to  transmit  forces  most  efficiently.  Since  the 
number  of  possible  material  choices  and  fiber  directions  is 
very  large,  suitable  methods  of  analysis  are  being  developed 
to  advance  the  methodology  for  minimizing  weight  in  the  de¬ 
sign  of  attachments  or  joints.  We  have  made  some  strides  in 
theoretical  analysis  and  intend  to  build  and  test  some  spe¬ 
cially  designed  specimens.  Bolted  attachments  involving 
only  a  single  bolt  are  being  looked  at  first. 

Our  theoretical  analysis  of  single  bolted  connections 
will  consider,  at  first,  not-yet-optimized  material  arrange- 
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merits  near  the  joint  consisting  of  ordinary  orthotropic 
graphi-te  epoxy  materials  and  plates  with ‘empty  holes.  This 
is  being  done  to  check  the  method  and  to  compare  our  efforts 
to  obtain  results  and  the  results  themselves  with  those  de¬ 
veloped  by  others.  The  next  case  considered  will  be  the 
loaded  hole  as  found  in  connections  such  as  the  attachments 
of  moving  control  surfaces.  The  manufacturing  problem  will 
be  kept  in  mind,  but  the  mechanical  efficiency  has  priority, 
at  least  for  these  first  stages. 

As  reported  in  previous  submissions,  the  Muser  design 
of  the  actuator  rib  of  the  B-727  elevator  was  based  on  the 
idea  of  using  fibers  running  largely  tangent  to  the  circum¬ 
ference  of  the  hole  through  which  the  forces  from  the  actuator 
rod  would  be  transmitted  into  the  elevator  structure.  It  is 
hoped  that  our  theoretical  work  v/ill  lead  to  a  new  design 
for  laminae  in  the  region  of  the  connector.  We  anticipate  a 
lay-up  pattern  with  fibers  tangential  to  the  circumference 
near  the  hole.  From  experience  and  from  published  reports 
(Ref.  12,  19,  20,  21  in  Table  V)  we  believe  that  this  arrange¬ 
ment  should  result  in  lower  stress  concentrations.  Our  work 
plan  has  the  objective  of  confirming  this  theoretically  and 
experimentally . 

It  should  be  noted  that  the  method  of  calculation  being 
developed  will  be  suitable  only  to  determine  stress  concen¬ 
tration  factors  in  the  linearly  elastic  range.  Since  bearing 
failure  occurs  outside  this  range,  calculations  valid  beyond 
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list  of  references  pertaining  to  bearing  strength  of  composites 


Frocht,  M.  M.  and  Ho  No  Hill,  "Stress  Concentration  Factors 
around  a  Central  Circular  Hole  in  a  Plate  Loaded  through  a 
Pin  in  the  Hole",  J.  Applied  Mechanics,  March  1940,  pp.  A5- 
A9o 

2.  Dalas,  Richard  N. ,  "Mechanical  Joints  in  Structural  Compos¬ 
ites",  Advances  in  Structural  Composites,  Vol.  12,  Society 
of  Aerospace  Material  and  Process  Engineers,  12th  National 
Symposium  and  Exhibit,  October  1967. 

3.  Dalas,  Richard  N.,  "Methods  of  Joining  Advanced  Fibrous 
Composites",  Composite  Materials;  Testing  and  Design,  ASTM 
STP  460,  1969,  pp.  381-392. 

4.  Schjelderup,  H.  C.  and  B.  H.  Jones,  "Practical  Influence  of 
Fibrous  Reinforced  Composites  in  Aircraft  Structural  Design", 
ASTM  STP  460,  1969,  pp.  285-306. 

5.  Waszczak,  J.  P.  and  T.  A.  Cruse,  "Failure  Mode  and  Strength 
Prediction  of  Anisotropic  Bolt  Bearing  Specimens",  J.  Com¬ 
posite  Materials,  Vol.  5,  1971. 

6.  Padawer,  G.  E.,  "The  Strength  of  Bolted  Connections  in 
Graphite/Epoxy  Composites  Reinforced  by  Colaminated  Boron 
Film",  ASTM  STP  497,  1972,  pp.  396-414. 

7.  Zurko,  M.  J.,  "Preliminary  Investigation  for  Developing  Op¬ 
timum  Fasteners  for  Joining  Composite  Materials",  Naval  Air 
Development  Center,  AD908930,  February  1973. 

8.  Greszczuk,  L.  B. ,  "Compressive  Strength  and  Failure  Modes 
of  Unidirectional  Composites",  ASTM  STP  521,  1973,  pp.  192- 
217. 

9.  Kim,  R.  Y.,  "Effect  of  Temperature  and  Moisture  on  Pin-Bear¬ 
ing  Strength  of  Composite  Laminates",  J.  Composite  Materials, 
Vol.  10,  April  1976,  pp.  149-155. 

10.  Hart-Smith,  L.  J. ,  "Bolted  Joints  in  Graphite-Epoxy  Compos¬ 
ites",  NASA  CR-144899,  June  1976. 

11.  Quinn,  W.  J.  and  F.  L.  Matthews,  "TheEffect  of  Stacking 
Sequence  on  the  Pin-Bearing  Strength  in  Glass  Fibre  Rein¬ 
forced  Plastic",  J.  Composite  Materials,  Vol.  11,  April 
1977,  139-145. 
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12.  de  Jong,  T. ,  "Stresses  around  a  Pin-loaded  Hole  in  an 
Elastically  Orthotropic  or  Isotropic  Plate",  J.  Compos¬ 
ite  Materials,  Vol.  11,  July  1977,  pp.  313-331. 

13.  Kubo,  J.  T. ,  "Joints",  Advanced  Composites  Design  Guide, 
Vol.  II,  Wright -Patter son  Air  Force  Base,  1977. 

14.  Collings,  T.  A.,  "The  Strength  of  Bolted  Joints  in 
Multi-Directional  CFRP  Laminates",  Her  Majesty's  Sta¬ 
tionary  Office,  London,  1977. 

15.  Wilkins,  D.  J. ,  "Environmental  Sensitivity  Tests  of 
Graphite-Epoxy  Bolt  Bearing  Properties",  ASTM  STP  617, 
1977,  pp.  497-513. 

16.  Eshwar,  V.  A.,  "Analysis  of  Clearance  Fit  Pin  Joints", 

J.  Mechanical  Science,  Vol.  20,  Pergamon  Press,  1978, 
pp.  477-484. 

17.  Fukuda,  H. ,  "Fastening  of  FRP",  Unpublished  Report, 
Aerospace  Research  Laboratory,  Tokyo  University,  1979. 

18.  Collings,  T.  A.,  "The  Use  of  Bolted  Connections  as  a 
Means  of  Joining  Carbon  Fibre  Reinforced  Plastics",  De¬ 
sign  with  Fibre  Reinforced  Materials,  Institute  of 
Mechanical  Engineers,  Conference  Publication,  1977-1979, 
pp.  29-32. 

19.  Puppo ,  A.  and  J.  Haener,  "Application  of  Micromechanics 
to  Joints  and  Cutouts",  U.  S.  Army  Aviation  Material 
Laboratories,  TR  69-25,  April  1969,  AD  688168. 

20.  Mueller,  J. ,  "Spannungsverteilung  an  gelochten  ueber 
einen  Bolzen  belasteten  Staeben  aus  kohlenstof faserver- 
staerktem  Kunststoff",  Swiss  Federal  Institut  of  Technol¬ 
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the  linear  region  would  probably  be  neceseary ,  and  suitable 
failure  criteria  would  have  to  be  established  before  there 
can  be  hope  for  good  agreement  between  analysis  and  test  re- 
suits.  We  expect  this  sort  of  work  to  wait  till  late  spring 
because  the  construction  of  the  test  specimens,  their  test¬ 
ing  and  the  linear  analytical  work  will  require  all  the  time 
available  until  then. 

4.  Supporting  Development  of  Mechanical  Joints 

The  work  of  graduate  student  Wonsub  Kim  and,  until 
August,  of  Research  Associate  Kiyoshi  Kenmochi,  in  this  pe¬ 
riod,  concentrated  on  the  interpretation  of  bearing  failure 
test  and  photoelastic  results  in  pin— loaded  composite  tension 
specimens.  This  followed  from  our  conclusion  that  designing 
a  mechanical  joint  to  fail  in  bearing  makes  the  best  use  of 
strength.  Reducing  fastener  diameter  increases  a  given 
joint's  tensile  strength  and  to  a  first  approximation  leaves 
shear  strength  unaffected,  but  bearing  stress  is  ultimately 
limiting  in  this  process.  For  this  reason,  we  have  concen¬ 
trated  on  bearing  strength  and  bearing  failure  mechanisms  in 
this  research. 

To  place  the  results  we  obtained  in  perspective,  a  lit¬ 
erature  survey  was  made  of  experimental  and  theoretical  v/ork 
published  in  this  area.  This  survey  has  shown  that  v/ithin 
the  substantial  literature  of  composite  materials,  relatively 
few  research  results  connected  with  m.echanical  joining  prob¬ 
lems  have  been  published.  Furthermore,  among  those  reported 
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works  which  do  exist  many  are  inconclusive  and  others  are 
not  mutually  consistent.  A  total  of  21  references  which 
seem  pertinent  are  listed  in  Table  V  on  pages  30-31  and  a  sum¬ 
mary  of  the  most  important  results,  from  the  point  of  view 
of  this  study,  follows.  ' 

Summary  of  the  Literature  Survey 

A.  General 

1)  There  are  relatively  few  publications 
available  dealing  with  bearing  failure 
in  composite  materials. 

2)  Among  such  publications  as  do  exist  the 
work  reported  is  roughly  90%  experimen¬ 
tal,  10%  theoretical. 

B.  Theoretical  V7ork 

1)  The  most  pertinent  work  appears  to  be 

* 

that  of  T.  de  Jong  (12)  and  V.  A. 

Eshwar  (16) .  Both  are  two  dimensional 
theories  and  assume  infinitely  rigid 
pins.  Neither  includes  the  effect  of 
friction,  comparisons  with  experiment 
nor  satisfies  boundary  conditions  ex¬ 
actly. 

2)  T.  de  Jong  shov/s  marked  effects  of  ply 
orientations  on  stress  distributions, 
but  provides  no  estimates  of  the  in¬ 
accuracies  introduced  by  failure  to 
meet  displacement  boundary  conditions 
on  the  hole. 


* 


Numbers  in  parentheses  in  this  section  correspond  to  refer¬ 
ence  numbers  in  Table  V,  pages  30-31. 
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3)  V.  S.  Eshwar's  work  is- limited  to  iso¬ 
tropic  materials  but,  by  using  an  inverse 
procedure  specifying  displacements  and 
then  solving  for  the  corresponding  loads 
on  the  hole,  shows  that  the  boundary  con¬ 
ditions  can  be  met  to  within  satisfactory 
accuracy.  Among  his  results  are  the 
following: 

a)  The  area  of  contact  between  pin  and 
plate  has  an  asymptotic  relation  to 
load  limited  by  Poisson's  ratio, 
for  example,  to  166°  (compared  to 
180°)  for  V  =  .3. 

b)  Non-linear  behavior  of  maximum 
stresses  with  load  are  exhibited 
(concave  downward)  as  a  result  of 
increasing  contact  area  with  load. 

c)  The  distribution  of  radial  stress, 
however,  changes  only  slightly 
with  load  and  is  well  approximated 

,  where 

9  is  the  angular  position  of  the 
c 

edge  of  contact  measured  from  the 
axis  of  symmetry. 

d)  The  effect  of  the  third  direction 
is  likely  to  be  great  near  the  hole 
and  suggests  caution  regarding  use 
of  2-D  theory  in  formulating  fail¬ 
ure  criteria. 

C.  Photoelastic  Studies 

Stress  concentration  factors  determined  by  Frocht 
and  Hill  (1)  were  found  .to  vary  as  follows: 


by 


a  =  a 
r  r 


max 
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1)  If  the  ratio  of  the  diameter  of  -the  hole 
to  the  width  of  the  plate,  D/t  increases, 
the  stress  concentration  factor  decreases. 

2)  If  the  clearance  between  pin  and  hole  in¬ 
creases,  the  stress  concentration  factor 
decreases . 

3)  If  the  number  of  pins  increases,  the 
stress  concentration  factor  decreases. 

4)  The  effect  of  lubrication  between  pin 
and  hole  is  measurable  but  very  small 
(thus  validating  the  assumption  of  no 
friction  between  pin  and  hole  in  mathe¬ 
matical  analyses) . 

Experimental  Work 

1)  All  referenced  v;orks  investigate  the 
effects  of  one  or  more  of  the  follow¬ 
ing  six  factors. 

a)  Fiber  orientation 

b)  Stacking  sequence 

c)  Ratio  of  hole  diameter  to  plate 
thickness  (D/t) 

d)  Clamping  force 

e)  Environmental  condition 

f)  Local  reinforcement  of  the  hole 

2)  All  tests  used  virtually  the  same  method, 
i.e.,  apply  load  to  the  pin  until  failure 
occurs. 

3)  Of  ten  references  reporting  the  effect  of 
fiber  orientation,  four  (3,  8,  15,  18)  re 
ported  the  effect  to  be  great,  five  (2,  9 
10,  14,  17)  pronounced  the  effect  small 
and  one  (13)  v/as  noncommital. 
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4)  Two  (9,  13)  of  three  investigations  reported 
environmental  factors  as  important  influ- 
enceSf  but  one  (15)  found  them  of  no  signifi¬ 
cance  = 

5)  All  three  of  those  reporting  stacking  se¬ 
quence  effects  (11,  14,  18)  concluded  it  is 
relatively  large  (But  two  of  these  were 
written  by  the  same  author.) . 

6)  Four  (10,  13,  14,  18)  out  of  six  papers  re¬ 
ported  clamping  force  as  having  a  large  in¬ 
fluence  on  bearing  strength.  The  remaining 
two  authors  (3,  7)  pointed  out  that  the 
significant  fact  is  the  existence  of  clamp¬ 
ing  force,  not  the  amount  of  such  force. 

7)  Five  investigators  tested  the  effect  of 
D/t;  one  (14)  reported  it  as  big,  two  (15, 

17)  as  small,  and  two  (2,  13)  found  their 
tests  inconclusive. 

8)  Three  authors  (2,  6,  10)  investigated  the 
influence  of  local  reinforcing  around  the 
hole  using  colaminated  third  materials, 

and  they  all  reported  the  effect  to  be  favor¬ 
able  and  very  large. 

9)  Tv;o  reports  (5,  13)  mention  the  effect  of 
clearance  between  pin  and  hole,  but  it  was 
only  calculated  using  a  Finite  Element 
Method . 

The  results  of  RPI's  testing  is  presented  in  this  sec¬ 
tion  in  some  detail  and  is  also  related  to  the  corresponding 
data  found  in  the  literature  search.  Table  VI  summarizes 
RPI  bearing  test  results  to  date.  Analysis  of  the  photo¬ 
elastic  test  results  hasn't  been  completed.  Some  initial 
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TABLE -VI 

RESULTS  OF  RPI  BEARING  TESTS 


★ 

Specimen 

Designation 

Specimen 
Thickness 
( inch) 

Load 

Maxim\m 

(lbs) 

Bearing 

Strength 

(psi) 

Strength  Ratio 
to  that  of 
Specimen  E 

** 

Type  A 

-1 

0.085 

1,240^^^ 

58,353 

0.82 

-2 

0.083 

1,500 

73,764 

0-94 

Type  B 

-1 

0.094 

1,610 

68,510 

0.97 

-2 

0.094 

1,725 

74,902 

0.95 

Type  C 

-1 

0.095 

1,445 

60,842 

0.86 

-2 

0.095 

1,625 

69,817 

0.89 

-3 

0.095 

1,120 

47,158 

0.60^^^ 

Type  D 

0.085 

1,400 

65,882 

0.93 

-2 

0.086 

1,310^^^ 

62,173 

0.79 

Type  E 

-1 

0.083 

1,470 

70,843 

1.00 

-2 

0.083 

1,600 

78,682 

1.00 

Photoelastic 

0.089 

2,750 

61,798 

0.72 

11  It 

0.095 

3,350 

70,526 

0.82 

II  II 

0.089 

3,100 

69,663 

0.81 

11  II 

0.085 

3,650 

85,882 

1.00 

*Specimen  Types  A-E  have  holes  with  h  inch  diameter. 

The  hole  diameter 

in  the  photoelastic  specimens  is  h  inch. 

'ft  ^ 

The  specimen  designations  are  interpreted  as  follows: 

Type  A 

{0°/90'’}8s 

Type  D 

{±45°/0''/90' 

»}4s 

Type  B 

{±45°/0°}6s 

Type  E 

{±45°/±22.5‘ 

V±67.5°/0°/90°}2s 

Type  C 

{±45°/90'’}6s 

drilled  with  high  speed  steel  drill  bit. 

(2) 

Found  to  be  premature  failure. 

^^^Ratioed  to  Case  E-2. 
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conclusions,  however,  can  be  drawn;  namely, 

A.  The  material  behaves  as  if  almost  perfectly  elastic 
until  just  in  the  vicinity  of  the  maximum  load. 
Through  both  the  unloading  and  loading  process  the 
charting  pen  follows  exactly  the  same  curve  of  load 
vs  displacement.  Even  after  the  unloading  and  load¬ 
ing  process,  the  photoelastic  fringes  were  identical 
to  those  observed  before  this  process.  The  behavior 
close  to  the  failure  load  could  not  be  checked  by 
these  means. 

B.  Stress  distributions  are  quite  different  from  one 
fiber  orientation  type  to  another  (see  Figure  10) . 

On  the  other  hand,  the  values  in  Table  VI  do  not 
always  confirm  the  photoelastic  insights.  "Scatter 
between  the  duplicate  specimen  pairs  is  large  enough 
to  almost  obscure  differences,  but  specimen  Type  E, 
which  is  closest  to  isotropy,  showed  the  highest 
strength;  Type  B  looks  next  best;  and  Types  A,  C  and 
£)  seem  roughly  the  same  and  least  good.  Since  hole 
diameter  to  specimen  width  varied  slightly  between 
the  photoelastic  and  other  specimens,  it  may  be  best 
to  consider  bearing  strength  averages  based  only  on 
non-photoelastic  tests.  Then  the  strength  of  lay-up 
Type  E  is  around  75  Ksi.  The  lowest  average  would 
be  that  of  Type  C  at  65  Ksi.  This  way  of  looking  at 
the  data  makes  the  maximum  variation  of  bearing 
strength  less  than  12%.  Hence,  contrary  to  what 
would  be  expected  from  photoelastic  results,  such  as 
shown  in  Figure  10,  the  pin-bearing  strength  does 
not  change  greatly  from  one  fiber  orientation  to 
another.  Tv/o  exceptions  to  this  statement  are  speci 
mens  which  failed  prematurely. 


Isochroma tic  fringe  patterns  for  various  fiber  orientations 
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Tables  VII  and  VIII  relate  the  RPI  results  to  those  in 
the  literature.  Table  IX  represents  an  attempt  to  draw  con¬ 
clusions  as  to  fiber  orientations  which  constitute  good  de¬ 
sign  for  bearing  strength  by  making  comparisions  selectively 
among  the  results  reported  in  the  literature  and  those  of 
the  RPI  tests  reported  here.  Since  materials  varied  from 
test  to  test;  bearing  failure  loads  are  normalized  to  par¬ 
ticular  ply  orientations.  VJhile  hardly  conclusive;  the  fol¬ 
lowing  tentative  conclusions  can  be  drawn: 

A.  Combinations  including  0°  and  ±45°  plies  have 
superior  bearing  strength.  (Comparisons  2,  3 
and  4  in  Table  V) 

B)  Having  90°  plies  provides  questionable  benefits. 
(Comparison  2  and  4) 

C)  Clamping  does  not  alter  the  relative  strength 
between  0/90  and  0/±45  plies.  (Comparison  2) 

D)  Clamping  does  seem  to  alter  the  relative 
strength  between  0/±45  and  90/±45  plies.  (Com¬ 
parison  4) 

E)  Changing  fiber  material  may  alter  the  relative 
strength  between  0/90  and  ±45  plies.  (Compari¬ 
son  1,  assuming  conclusion  C)  holds  also  for 
those  ply  comparisons) 

An  interesting  conclusion  from  our  testing  is  that  the 
bearing  strength  does  not  seem  to  be  a  strong  function  of  the 
ratio  of  the  diameter  of  the  hole  to  the  thickness  of  the 
plate.  This  is  contrary  to  results  presented  in  Reference 
14.  This  D/t  ratio  for  photoelastic  specimens  was  roughly 
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BEARING  STRENGTH  WITH  VARIOUS  FIBER  ORIENTATIONS 
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TABLE  TX 


NORMALIZED  BEARING  STRENGTH  COMPARISONS  AMONG  REFERENCE  DATA 
—  AND  FROM  ~~ 

RPI  WORK  TO  DETERIvIINE  THE  EFFECT  OF  FIBER  ORIENTATIONS 


Comparison 

Number 


Fiber 

Orientation 


1. 

* 

Ref.  No.  —  14 

17 

Mat'l.  —  G/E 

Glass/E 

{0/90}s 

1.0 

1.0 

{±45}s 

1.04 

1.2 

2. 

* 

Ref.  No.  —  14 

RPI 

Mat '  1 .  —  G/E 

G/E 

{0/90}s 

1.0 

1.0 

{0/±45}s 

1.12 

1.17 

3. 

* 

Ref.  No.  —  14 

** 

6 

Mat'l.  —  G/E 

G/E 

{±/45}s 

1.0 

1.0 

{0/±45}s 

1.08 

1.7 

4. 

Ref .  No . 

Mat ' 1 . 

9 

—  G/E 

*** 

15 

G/E 

RPI 

G/E 

{0/±45}s 

1.0 

1.0 

1.0 

{90/±45}s 

0.88 

0.51 

0.89 

{0/90/±45}s 

1.06 

— 

0.97 

*3.2  Ksi  clamping;  **200  lb  clamping;  "5  lb-in  on  a  screv/" 
clamping. 


5.6  and  for  the  bearing  coupons  was  2.8,  yet  the  correspond¬ 
ing  bearing  strengths  were  quite  close.  With  such  ratios, 
according  to  Reference  14,  the  bearing  strength  of  the  photo¬ 
elastic  specimen  should  be  lower  than  the  bearing  test  coupon 
by  at  least  50%,  but  this  series  of  tests  at  RPI  shows  that 
the  photoelastic  specimens  have  strengths  which  are  even 
slightly  higher. 

5.  Failure  Mechanisms 

Some  of  the  specimens  tested  to  failure  in  our  program 
were  cut  along  the  center  line  passing  through  the  two  holes 
so  as  to  allow  microscopic  investigation  through  the  speci¬ 
men  thickness  in  the  area  of  the  failure. 

After  very  careful  polishing,  the  cut  edges  were  in¬ 
spected  both  on  the  failed  area  of  the  hole  as  well  as  on 
the  unfailed  side.  Microscope  magnification  factors  of  25 
and  50  were  used.  With  the  magnification  factor  of  50,  the 
cracks  on  the  damaged  area  as  well  as  each  layer  of  fiber 
orientation  could  be  seen  clearly.  Figure  11  shows  the 
cracks  on  failed  specimens  with  several  typical  fiber  orien¬ 
tations.  These  pictures  were  taken  with  a  magnification 
factor  of  25. 

It  seems  clear  that  the  cracks  started  at  an  arbitrary 
point  on  the  loaded  edge  of  the  hole  and  ran  across  the  lay¬ 
ers,  even  through  the  0®  layers,  in  a  direction  roughly  ±45° 
to  the  laminae.  This  kind  of , crack  has  been  observed  often 


Cracks  within  the  laminates  following  bearing  failure  (Magnif ication25) 
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in  tests  of  brittle  isotropic  materials  loaded  in  compres¬ 
sion.  If  these  cracks  are  along  shear  lines  as  they  appear 
to  be  and  as  they  are  in  compression  tests  of  isotropic  ma¬ 
terials,  then  shear  in  a  direction  transverse  to  the  plane 
of  the  plate  should  be  investigated  and  regarded  as  an  im¬ 
portant  factor  in  the  bearing  failure  of  composite  materials. 

6.  Special  Delamination  Studies 

The  effect  of  delamination  on  the  bearing  strength  of 
a  composite  material  seems  to  be  much  more  severe  than  it 
is  on  other  material  allowables.  In  a  brief  investigation 
of  this  effect  two  additional  specimens  were  prepared  using 
exactly  the  same  procedures  as  those  used  in  all  the  other 
tests.  In  this  case,  hov/ever,  the  holes  were  purposefully 
drilled  with  a  high  speed  steel  drill  bit  instead  of  a  car¬ 
bide  tipped  drill  bit.  Every  dimension  and  process  of  prep¬ 
aration  was  identical  except  the  drill  bit  used.  As  expected 
these  holes  were  flawed  to  a  greater  extent  than  those  drill¬ 
ed  with  a  carbide  tipped  bit.  The  thickness  of  the  plates 
around  these  holes  was  increased  by  up  to  one  thousandth  of 
an  inch.  This  indicates  that  delaminations  were  induced  by 
drilling  with  the  high  speed  steel  bit.  These  holes  were  in¬ 
vestigated  further  using  ultrasonic  NDE  techniques. 

Figure  12,  showing  ultrasonic  transmission  loss  images 
of  the  delamination  around  holes  drilled  vfith  a  high  speed 
steel  bit  and  with  a  carbide  tipp®^  bit,  shov7s  the  former  to 
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be  more  extensive.  Note  that  the  technique  does  not  yet 
allow  interpretation  of  delamination  as  a  function  of  depth, 
but  only  of  area  around  the  hole  in  the  plane  of  the  lami“ 
nae. 

These  two  specimens  were  tested  using  the  same  process 

as  described  above.  A  mishap  failed  one  specimen  of  lay-up 

Type  A  ({0°/90°},)  )  before  useful  data  could  be  taken.  Spec 

s 

imen  Type  C  ({±45?/90°}g^) ,  however,  was  tested  successfully 
It  shov/ed  much  less  strength  than  the  same  specimen  with 
higher  hole  quality,  as  expected.  This  result  is  also  in¬ 
cluded  in  Table  II. 

These  special  tests  lead  to  the  following  conclusions: 

1)  Holes  drilled  with  carbide  tipped  bits 
have  delamination  like  that  caused  by 
high  speed  steel  bits,  only  less  of  it. 

2)  The  amount  of  delamination  caused  by 
drilling  affects  the  bearing  strength 
of  a  joint  substantially. 

3)  Preventing  or  compensating  for  delami¬ 
nation  in  some  way  can  improve  the  bear¬ 
ing  strength  of  a  joint  greatly.  The 
large  and  favorable  effect  of  clamping 
seems  to  stem  from  this  phenomenon. 


PART  II 


CAPGLIDE  (Composite  Aircraft  Program  Glider) 
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CAPGLIDE  (Composite  Aircraft  Program  Glider) 

(R.  J.  Diefendorfy  H.  J,  Hagerup,  G.  Helwig) 

CAPGLIDE  is  an  undergraduate  program  to  design,  build 
and  test  advanced  composite  structures.  Students  will  ob™ 
direct  "hands~on"  experience  in  advanced  composite 
structures  which  can  serve  as  a  springboard  for  the  more 
sophisticated  CAPCOMP  projects.  In  dealing  with  the  design 
of  a  complete  vehicle,  the  effect  of  any  given  change  on 
other  aspects  must  be  dealt  with.  In  this  way  the  project 
g^^go  nsguires  students  majoring  in  aeronautical,  mechanical 
and  materials  engineering  to  interact  in  much  the  same  way 
as  they  do  in  industry. 

An  ultra-light  sailplane  was  selected  as  the  first 
demonstration  project  because  a  full  scale  flight  vehicle 
would  maximize  student  interest  and  would  be  of  relative 
simplicity  and  lov7  cost  to  build.  A  conventional  layout 
monoplane  with  three-axis  control  resulted  in  the  following 
estimated  performance: 

1)  Stall  speed,  15  knots, 

2)  Best  glide  ratio,  17  and 

3)  Minimum  sink  rate,  2.0  feet  per  second. 

While  the  glide  ratio  of  the  ultra-light  sailplane  is  simi¬ 
lar  to  that  of  post  VJorld  War  II  utility  gliders,  the  more 
important  sink  rate  is  in  the  range  of  standard  class  sail¬ 


planes  . 


Wet-Lay-Up  Mold  Technique 
Figure  13 
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Before  the  cap-strips  of  the  wing  spars  were  manufac¬ 
tured,  a  test  spar,  which  was  built  earlier,  was  loaded  to 
destruction.  A  test  loading  jig  was  hung  from  a  crane  and  was 
used  to  introduce  the  forces  at  four  points  into  the  canti¬ 
levered  spar  (see  Figure  14) . 


Figure  14  Spar  Test 

The  cantilever  point  was  achieved  by  clamping  with  the  con¬ 
nection  from  the  other  wing.  The  wing -connection  was,  there¬ 
fore,  tested  in  the  same  process.  Special  guides  were  neces¬ 
sary  to  prevent  the  I-beam  spar  from  twisting,  and  these  can 
also  be  seen  in  Figure  14.  In  the  first  attempt,  the  loading 
jig  failed  at  50%  of  the  ultimate  load  (i.e.,  design  load 
with  safety  factor  of  1.5),  and  the  test  spar  was  damaged  and 
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had  to  be  repaired.  The  next  test  was  successful.  The  spar 
failed  in  the  connection  with  a  safety  factor  slightly  over 
1.5.  The  nature  of  the  failure  was  graphite  capstrips  split¬ 
ting  under  shear  forces.  The  actual  wing  connection  will  be 
reinforced  to  prevent  shearing-off  of  the  I-beam  flanges.  A 
still  higher  safety  factor  can,  therefore,  be  expected.  A 
3-foot  wing  section,  previously  tested  for  wing-skin  buckling, 
was  tested  again  at  an  elevated  temperature  of  about  125°F. 
Significant  stiffness  loss  was  obvious.  Accordingly,  the 
wing  was  reinforced  with  two  additional  webs  from  the  root 
to  the  inboard  end  of  the  ailerons,  which  added  some  extra 
weight. 

The  tail  was  built  with  single-sided,  Kevlar-stabilized 
foam  sheets.  These  sheets  were  then  bent  to  the  rib  frame¬ 
work.  This  technique,  again,  was  simple  and  effective  in 
producing  true— to— tolerance  airfoil  shapes  with  an  extremely 
lightweight  structure. 

From  the  following  Figures,  Id— 18,  the  construction 
details  and  materials  actually  used  in  fabricating  the  glider 
can  be  seen.  The  actual  weight  of  the  parts  is  also  listed. 
Weights  of  major  components  sub-assemblies  are  shown  in  Table 
X.  Figures  19-21  are  photographs  which  show  various  assem¬ 
blies  of  the  first  CAPGLIDE  aircraft,  RP-1,  in  different 
stages  of  completion;  Figure  19  is  the  tail  assembly,  Figure 
20  is  a  semi-span  wing  substructure  and  Figure  21  is  the 
fuselage. 
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_ Item _ 

Span 

Length 

Height 

Wing  area 

Aspect  ratio 

Wing  section  profile 

Vertical  stabilizer 
Area 
Span 
Profile 

Horizontal  stabilizer 
Area 
Span 
Profile 

Performance  with  120  kg 
Stall  speed 
Minimum  sink  rate 
at 

Maximum  glide-ratio 

at 

Empty  weight 

Tail  control  effectivenesses 
V/V*  Horizontal 
V/V*  Vertical 


Units 

m  (f  t) 

11.5 

(37.7) 

m  (f  t) 

5.78 

(18.9) 

m  (ft) 

1.28 

(4.2) 

2  2 
m  (ft  ) 

12 

11 

(129.0) 

FX-63-137 

2  2 
m  (ft  ) 

0.72 

(7.8) 

in(ft) 

1.2 

FX-L-111-142 

(3.9) 

2  2 
m  (ft  ) 

1.8 

(19.4) 

m  (ft) 

3.0 

FX-L-111-142 

(9.8) 

km/h (mph) 

36 

(22.5) 

m/s (ft/s) 

0.62 

(2.0) 

km/h (mph) 

40 

20 

(25.0) 

km/h (mph) 

52 

(32.5) 

kg (lbs) 

52.9 

(116.3) 

0.41 

0.057 


Figure  16  RP-1 
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Part 

Materials 

Weight  (g) 

2  Spar  capstrips 

Graphite 

1,600 

Main  web 

Kevlar,  foam 

700 

17  Ribs 

Glass,  foam 

.1,020 

2  Skins 

Kevlar,  foam 

6,580 

Front  web 

Glass,  foam 

14  0 

Reinforcement  webs 

Glass,  foam 

340 

Back  web 

Glass,  foam 

83 

Root  rib 

Graphite,  foam,  spruce 

750 

Wing  connection 

Graphite,  plywood, 

foam 

500 

Aileron  tube 

Graphite 

220 

Hinges 

Graphite 

100 

Wing  tips 

Kevlar,  foam 

900 

Leading  edge 

Resin  for  bonding 

Balsa 

Total 

500 

13,433 

2,600 

=■  16,033 

Resin  for  bonding 
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RENSSELAER  POLYTECHNIC 
INSTITUTE  NOVEMBER  1979 


Part 

Materials 

Weight  (g) 

4  Boom  capstrips 

Graphite 

3,000 

4  Webs 

Kevlar,  foam 

3,400 

Stiffeners 

Kevlar,  foam 

200 

Fuselage 

Kevlar 

900 

Pilot  beam 

Graphite,  balsa 

1,300 

Landing  skid 

Glass,  spruce 

1,100 

4  Wing  supports 

Spruce 

320 

Push  Rods 

Aluminum 

500 

Aileron  Tube 

Graphite 

80 

Cranks 

Graphite 

100 

Control  stick 

Graphite 

120 

4  Landing  struts 

Graphite,  balsa 

600 

2  Connections  blocks 

Spruce 

140 

Torsion  box 

Graphite,  foam 

810 

3  Boom  connections 

Graphite,  balsa 

200 

8  Filler  blocks 

Spruce,  balsa 

280 

13,050 

Resin  for  bonding 

2,500 

Total 

=  15,550 
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Figure  18 

TAIL-CAPGLIDE 

RENSSELAER  POLYTECHNIC 
INSTITUTE  NOVEMBER  197® 


Part 

Materials 

Weight  (g) 

Skins 

Kevlar,  foam 

3,000 

24  Ribs 

Foam 

144 

Web  (horizontal) 

Kevlar,  foam 

192 

Web  (vertical) 

Kevlar,  foam 

96 

2  Capstrips 

Graphite 

225 

4  Rudder  Webs 

Foam 

54 

11  Hinges 

Graphite 

220 

Rudder  connection 

Graphite,  tube 

150 

Leading  edge 

Balsa 

135 

Wing  tips 

Balsa 

90 

2  T-Webs 

Kevlar,  foam 

128 

Tail  connection 

Graphite 

50 

4,484 

Adhesive 

600 

Total 

=  5,084 
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TABLE  X 

CAPGLIDE :  RP-1  ACTUAL  TOTAL  WEIGHTS 


(kg) 

(lbs) 

Wings 

32o2 

o 

00 

Tail 

5.1 

11.2 

Fuselage 

15.6 

34.3 

Total  Weight 

52.9 

116.3 

Tail  Assembly 


Figure  19 


PART  III 


COMP AD  (Computer  Aided  Design) 
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COMPAD  (Computer  Aided  Design) 

(L.  J.  Feeser  and  M.  S.  Shephard) 

■  During  this  reporting  period,  emphasis  within  the  com¬ 
puter  aided  design  portion  of  the  composites  project  has 
been  put  on  enhancing  RPI's  finite  element  capabilities. 

This  work  is  primarily  concerned  with  three  activities: 

(1)  extension  and  maintenance  of  the  finite  element  analysis 
programs  available  on  RPI's  computer  systems,  (2)  initial 
development  of  a  highly  interactive  finite  element  prepro¬ 
cessor  and  (3)  linking  the  finite  element  analysis  program, 
SPAR,  to  an  array  processor. 

1.  Finite  Element  Analysis  Programs 

Work  has  begun  on  the  development  of  a  nev/  in-house 
finite  element  analysis  program.  This  work  is  being  carried 
out  primarily  as  a  series  of  term  projects  by  a  few  selected 
students  in  RPI's  finite  element  courses.  This  new  program 
is  being  developed  in  a  highly  modular  form  so  that  new  cap¬ 
abilities  can  be  added  with  a  minimum  of  effort. 

Several  other  finite  element  method  analysis  programs 
are  available  to  us;  in  fact,  enhancement  of  RPI's  finite 
element  analysis  capability  has  been  accomplished  mainly  by 
the  acquisition  of  additional  finite  element  analysis  codes. 
Table  XI  lists  the  finite  element  programs  currently  avail¬ 
able  along  with  a  brief  description  of  their  capabilities. 
Nevertheless,  we  anticipate ' that  at  some  point  in  the  future 
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TABLE  -  XI 

list  of  finite  element  codes  available  on  rpi's  computer  systems 


Program 


SPAR  PRIME 


SAP  IV 


NONSAP 


IBM, 

PRIME 

IBM 


ASHSD2  IBM 

STRUDL  IBM 

* 

ansys  prime 


Capability 


Static  and  dynamic  analysis 
of  both  linear  and  nonlinear 
structural  systems.  Especi¬ 
ally  well  suited  for  compos¬ 
ites  . 

Static  and  dynamic  analysis 
of  linear  three-dimensional 
structural  systems. 

Same  as  SAP  IV  but  with  non¬ 
linear  capa bil r t le s o  Nonlin 
earities  may  be  due  to  large 
displacements,  large  strains 
and  nonlinear  material  behav¬ 
ior  . 

Evaluates  the  time  dependent 
displacements  and  stresses  of 
complex  ax i symmetric  struc¬ 
tures  subjected  to  any  arbi¬ 
trary  static  or  dynamic  load¬ 
ing  or  base  accelerations. 

Static  and  dynamic  analysis 
of  linear  structural  system. 

A  wide  variety  of  capabili¬ 
ties  are  available  in  this 
program  including: 

a)  static  and  dynamic  analy¬ 
sis  of  linear  and  nonlin¬ 
ear  systems  (elastic, 
plastic,  creep,  swelling 
and  large  deflections 
handled) , 

b)  steady-state  or  transient 
heat  transfer  analysis 
(conduction,  convection 
and  radiation)  and 

c)  some  fluid  and  electrical 
capabilities . 


*By  contract  the  use  of  ANSYS  is  limited  to  "small"  problems 
for  course  work  and  student  projects.  It  will  not  be  avail¬ 
able  until  Fall  '79. 
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we  will  need  to  develop  some  new,  then  unavailable  capabil- 
ites,  and  it  will  be  far  more  efficient  to  add  these  capabil- 
ites  to  a  general  program  that  we  know  intimately,  than  it 
would  be  to  modify  one  of  the  other  existing  programs  or  to 

a  complete  program.  An  additional  enhancement  which 
will  be  discussed  in  conjunction  with  the  array  processor 
work  is  the  creation  of  a  version  of  SPAR  that  will  run  on  a 
PDP  11/40  minicomputer. 

2.  Preprocessor  Development 

The  finite  element  preprocessor  currently  under  develop¬ 
ment  will  handle  virtually  all  possible  two-dimensional  to¬ 
pologies.  Extensive  use  of  graphical  input  techniques  and  a 
unified  overall  design  will  result  in  a  preprocessor  that  is 
both  efficient  and  simple  to  use. 

The  development  of  a  finite  element  model  using  a  general 
preprocessing  capability  can  be  thought  of  as  a  three  stage 
process.  The  first  and  most  general  phase  is  defining  and 
discretizing  the  problem  geometry.  The  second  phase  is  de¬ 
fining  problem-specific  attributes  such  as  material  proper¬ 
ties,  boundary  conditions  and  loads.  The  third  phase  is 
linking  the  preprocessor  with  the  particular  analysis  pro¬ 
gram  selected  to  solve  the  model.  With  this  approach,  only 
one  general  program  for  defining  and  discretizing  problem 
geometry  is  required,  since  this  process  is  the  same  for  all 
finite  element  models.  This  program  can  then  be  linked  to 
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one  attribute  editor  for  each  class  of  problem  to  be  con¬ 
sidered.  For  example,  one  attribute  editor  may  be  used  for 
stress  analysis  problems,  another  for  heat  transfer  analysis 
and  a  third  for  fluid  flow  calculations.  The  last  stage, 
the  link  to  an  analysis  routine,  is  little  more  than  a  sim¬ 
ple  reformatting  program?  it  is  the  only  portion  of  the 
preprocessor  program  that  will  have  to  be  repeated  for  each 
new  analysis  routine  considered. 

A  two-dimensional  preprocessor  has  been  developed  to 
the  point  that  basic  problem  geometries  can  be  defined  and 
element  meshes  can  be  generated  using  a  set  of  mesh  gener¬ 
ating  capabilities.  The  operation  of  the  preprocessor  will 
be  demonstrated  in  terms  of  a  simple  example. 

The  first  step  is  the  definition  of  a  background  plan¬ 
ning  grid  which  can  be  thought  of  as  a  piece  of  electronic 
graph  paper.  This  planning  grid  is  used  later  to  obtain  nu¬ 
merically  exact  input  with  interactive  input  devices  by  "hit 
testing"  on  the  intersection  of  grid  lines.  In  addition  to 
defining  and  scaling  a  basic  rectangular  grid,  interactive 
features  allow  the  user  to  add  and  delete  additional  vertical 
and  horizontal  grid  lines  at  user-specified  levels  and  to 
move  the  origin  to  any  grid  intersection  to  produce  a  final 
planning  grid  such  as  the  one  shown  in  Figure  22. 

The  second  step  is  to  define  the  subregion  boundary 
curves  to  be  used  by  the  mesh  generators.  Curves  are  defined 
interactively  by  first  selecting  the  type  of  curve  to  be 
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generated  and  then  inputting  the  appropriate  number  of  "key 
points"  (either  by  use  of  the  light  pen,  thumb  wheels  or 
keyboard)  along  the  curve »  The  curve-generating  options 
include  single  points,  lines,  circles,  circular  arcs  and 
second  and  third  degree  polynomials  computed  by  Lagrangian 
interpolation.  A  specified  number  of  "key  nodes"  are  then 
placed  along  the  curve.  The  key  nodes  are  the  nodal  loca¬ 
tion  used  by  the  mesh  generator  to  define  the  resulting  mesh. 
Since  the  number  and  location  of  the  key  nodes  define  the  re¬ 
sulting  mesh,  options  are  available  to  specify  both  the  num¬ 
ber  of  key  nodes  to  be  placed  on  a  curve  and  also  a  weighting 
parameter  that  forces  the  nodal  spacings  to  vary  from  uniform 
along  the  length  of  the  curve  to  being  biased  toward  either 
end.  Figure  23  shows  the  interactive  input  of  a  circular 
arc  by  three  points  along  the  arc  while  Figure  24  shows  the 
final  set  of  subregion  boundary  curves  for  the  example  prob¬ 
lem. 

After  the  subregion  boundary  curves  have  been  defined, 
the  geometry  of  the  object  has  been  determined  and  the  element 
mesh  can  be  generated.  The  element  mesh  is  generated  one  sub- 
region  at  a  time  by  employing  one  of  several  available  mesh 
generators . 

Currently,  three  mesh  generators  have  been  implemented. 
The  generators  are  computationally  efficient  and  will  produce 
well-conditioned  grids  in  highly  distorted  regions.  Each 
generator  is  based  on  a  mathematical  mapping  that  defines 
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internal  node  points  in  terms  of  the  appropriate  boundary 
key  nodes  by  linear  blending  interpolants  in  one,  two  or 
three  directions*.  After  the  locations  of  the  internal  nodes 
have  been  defined,  elements  can  be  specified  by  appropriately 
connecting  the  nodes. 

The  mesh  generators  are  invoked  by  first  interactively 
indicating  the  subregion  boundary  curves  that  define  the 
sides  of  a  subregion,  selecting  an  element  type  and  then  in¬ 
dicating  the  mesh  generator  to  be  used.  Considering  the  left 
hand  subregion  of  the  example  problem,  the  user  indicates 
that  the  first  three  sides  are  one  curve  each,  while  the 
fourth  side  (the  bottom  of  the  subregion)  is  the  concatena¬ 
tion  of  two  curves.  The  user  then  indicates  that  he  wants  a 
four-sided  subregion  generated  by  the  linear  blending  func¬ 
tions  using  four  noded  quadrilateral  elements  to  define  the 
mesh.  The  resulting  mesh  is  shown  in  Figure  25  where  a  zoom 
feature  has  been  used  to  look  more  closely  at  that  subregion. 
Triangular  elements  can  also  be  generated  as  shown  in  Figure 
26.  Figure  27  shows  triangular  elements  generated  in  a 
three-sided  region. 

The  upper  right  hand  portion  can  be  generated  as  a  single 
subregion  by  mapping,  via  the  linear  blending  functions,  from 
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the  circle  defining  the  hole  to  the  three  curves  defining 
the  outside  of  the  subregion.  Figure  28  shov/s  a  four  element 
deep  mesh,  biased  towards  the  hole,  generated  in  this  manner. 
Figure  29  shov/s  the  completed  mesh  for  the  example  problem. 

3.  Array  Processor  for  SPAR 

As  discussed  in  the  last  semi-annual  progress  report, 
preliminary  feasibility  studies  have  indicated  the  potential 
benefits  of  using  an  array  processor  in  conjunction  with  the 
SPAR  finite  element  analysis  program.  Work  is  currently 
under  way  to  link  SPAR  with  a  Floating  Point  Systems  AP120 
array  processor.  Since  this  array  processor  is  operated 
from  a  POP  11/40  mini-computer,  effort  has  concentrated  on 
creating  a  modified  version  of  SPAR  to  run  on  the  11/40.  As 
this  work  is  being  conducted,  areas  of  code  where  an  array 
processor  may  be  used  are  being  identified  and  considered  as 
candidates  for  replacement  by  array  processor  operations. 
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PART  IV 

INSURE  (Innovative  and  Supporting  Research) 


Progress  on  composites  research  is  reported  in  the  indi¬ 
vidual  write-ups  on  the  succeeding  pages  in  the  following 
areas ; 

Advanced  Analysis  Methods  for  Composite  Structures 
Ultrasonic  Non-Destructive  Testing  Developments 
Physical  Properties  of  Epoxy  Resins  and  Composites 
Fatigue  in  Composite  Materials 

Transverse  Thermal  Expansion  of  Carbon/Epoxy  Composites 
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ADVANCED  ANALYSIS  METHODS  FOR  COMPOSITE  STRUCTURES 

Senior  Investigator:  E.  J.  Brunelle 


During  the  reporting  period  three  diverse  areas  of  re¬ 
search  were  undertaken.  The  first  two  areas  appear  to  be 
relatively  short  term  ventures  even  though  one  of  them  has 
become  diffuse  enough  to  lead  to  a  number  of  sub-areas  that 
are  each  producing  interesting  results.  The  third  area  is  a 
more  ambitious  undertaking  and  is  considered  to  be  a  longer 
term  effort  with  potentially  promising  and  useful  results. 
Interestingly  enough,  research  in  the  first  two  areas  was 
initiated  only  after  assuming  that  the  results  should  have 
been  in  the  literature;  however  no  trace  of  them  was  found, 
and  it  v/as,  therefore,  deemed  worthwhile  to  spend  some  time 
filling  in  the  theoretical  gaps. 

The  first  area  of  research  delineates  the  lav/s  governing 


the  variation  of  the  properties  of  an  orthotropic  lamina,  of 
an  anisotropic  lamina  and  of  a  laminate  comprised  of  ortho- 
tropic  laminas  with  changes  in  the  orientation  of  the  coordi- 
nate  axes.  A  nondimensional  number  B  =  completely  governs 

these  rotational  properties  of  an  orthotropic  lamina  and  de¬ 
pending  on  whether  B>1  or  B<1,  "regular"  or  "anomalous"  rota¬ 
tional  properties  are  observed.  This  is  not  simply  an  "aca¬ 
demic"  issue  since  typical  values  of  B  for  boron/epoxy, 
graphite/epoxy  and  glass/epoxy  are  .99,  1.01  and  1.3  respec¬ 
tively.  Besides  obtaining  detailed  rotational  information  on 
the  Q • j  (elements  of  the  reduced  stiffness  matrix  in  the 


transformed  stress-strain  relation)  as  influenced  by  B,  a 
simple  inequality  has  been  derived  that  identifies  a  new  com¬ 
posite  material  as  being  "regular"  or  "anomalous".  This  in¬ 
equality  (see  Figure  30)  should  be  of  particular  interest  to 
the  physical  chemists/materials  engineers  who  are  involved 
with  developing  new  fibrous  composite  materials. 

The  second  area  of  research  traces  its  origins  to  a  de¬ 
sire  to  reduce  the  number  of  elastic  constants  in  any  given 
problem  to  a  smaller  number  of  "cluster  constants"  which  can 
be  used  more  efficiently  in  parameter  studies  of  any  given 
problem.  For  example,  if  a  problem  involving  four  elastic 
constants  could  be  reduced  to  a  problem  involving  only  one 
"cluster  constant",  physical  insight  into  the  problem  becomes 
much  clearer,  more  general  solutions  are  usually  possible. 
Furthermore,  what  was  originally  an  enormous  computational 
effort  (Note  that  for  a  parametric  study  of  the  above  cited 
generic  problem,  sufficient  sets  and  combinations  of  the  four 
constants  must  be  considered.)  can  be  reduced  to  a  one  para¬ 
meter  family  of  solutions. 

"Cluster  constants"  have  been  derived  by  employing  af¬ 
fine  transformations  on  the  dependent  variables  as  well  as 
the  independent  variables.  This  method  has  been  quite  suc¬ 
cessful  and  has  produced  some  diverse  results  which  follow: 

1)  A  composite  shell  problem  that  was  posed,  but 
not  solved  correctly,  almost  two  decades  ago 
in  a  classic  paper  by  Pister,  Dong  and  Taylor 
is  now  easily  solved.' 
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2)  The  dependence  of  shear  deformation  effects 
on  boundary  conditions  has  been  shown  drama¬ 
tically  by  a  general  solution  to  a  composite 
plate  strip  problem»  The  result >  which  is 
well  documented  in  earlier  work  for  (simpler) 
transversely  isotropic  beams  and  plates,  is 
that  the  effects  of  shear  deformation  are 
exacerbated  by  any  form  of  increased  bound 
ary  restraint;  obviously  this  is  not  the 
trend  one  would  prefer  to  see. 

3)  The  specially  orthotropic  plate,  which  pro¬ 
vides  important  bounds  for  composite  plate 
problems  (upper  bounds  for  buckling,  etc.), 
can  be  reduced  to  a  one  parameter  family  of 
solutions.  Master  curves  can  thus  be  pre¬ 
sented  for  any  and  all  static,  dynamic  and 
instability  problems  associated  with  the 
specially  orthotropic  plate. 

The  third  area  of  research  seeks  a  systematic  method  of 
generating  hierachies  of  (elastic)  composite  beam,  plate  and 
shell  equations  that  govern  the  static,  dynamic  and  stability 
behavior  of  realistic  composite  materials  structures.  In 
this  development,  a  proper  accounting  is  to  be  made  of  (i) 
conservative  and  non— conservative  arbitrary  initial  stress 
states,  (ii)  the  accompanying  initial  displacement  gradient 
field  [this  has  been  almost  completely  neglected  in  the  lit¬ 
erature  of  even  single  layer  structures]  and  (iii)  arbitrary 
states  of  initial  [stress-free]  displacements  so  that  the 
effects  of  "imperfection  sensitivity"  are  an  integral  part 
of  the  formulation. 
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The  progenitor  of  this  method,  as  developed  by  the  pres¬ 
ent  investigator,  will  reproduce  all  the  standard  results  of 
structural  mechanics  (with  some  interesting  addtional  features) 
when  the  simplest  possible  displacement  assiamptions  of  the 
hierachy  are  introduced  into  a  non-classical  virtual  work 
theorem.  The  results  of  this  research  should  be  considered 
as  a  serious  candidate  to  provide  the  basic  equations  to  be 
used  in  the  next  generation  of  finite  element  programs  for 
routine  structural  analysis.  It  may  be  desireable  to  give 
these  finite  element  method  (FEM)  programs  self —diagnostic 
capability  such  that  the  programs  evaluate  the  input  and 
automatically  select  the  proper  hierachy  of  equations  that 
v/ill  most  properly  and  most  economically  solve  the  given 
structural  problem.  In  other  words,  the  self -diagnostic 
capability  should  be  designed  to  prevent  "overkill"  or 
"underkill"  in  the  solution  of  any  given  problem. 
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ULTRASONIC  NON -DESTRUCTIVE  TESTING  DEVELOPMENTS 

Senior  Investigators:  H.  F.  Tiersten 

P.  K.  Das 

Since  the  last  report  was  submitted,  trapped  energy  mode 
focusing  transducers  have  been  fabricated  by  plating  a  number 
of  concentric  rings  on  a  PZT-7A  plate.  The  concentric  ring 
structure  acts  as  a  Fresnel  lens  and  focuses  the  ultrasonic 
waves.  Two  transducers,  operating  in  the  frequency  range 
centered  at  2.5MHz,  have  been  fabricated.  The  first  order 
diffraction  pattern  from  a  laser  probe  is  being  used  to  study 
the  radiation  pattern  in  water  at  various  distances  from  the 
transducer.  From  the  near  field  radiation  pattern  it  has 
been  confirmed  that  the  acoustic  isolation  between  the  rings 
is  excellent. 

A  doubly  focused  acoustic  imaging  system,  which  incor¬ 
porates  focusing  transducers  at  both  the  transmitting  and 
receiving  ends,  is  being  used  for  imaging  of  composite  mate¬ 
rials.  Scanning  is  performed  using  a  Z-80  microprocessing 
system  which  also  stores  and  displays  the  image.  Images  of 
the  ultrasonic  transmission  through  composite  materials, 
which  were  prepared  for  tensile  testing,  are  being  obtained 
(See  Part  I,  CAPCOMP,  section  4).  Preliminary  results  com¬ 
pare  well  with  those  obtained  using  other  testing  procedures 
and  reveal  certain  features  which  have  never  been  observed  by 
means  of  non-destructive  testing  procedures. 
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The  analysis  of  mode  coupling  in  thickness-extensional 
trapped  energy  mode  transducers  with  strip  electrodes  has 
been  completed.  The  frequency  spectra  occur  in  closed  loops 
due  to  the  large  coupling  in  PZT-7A.  The  analysis  uncovered 
modes  with  numerous  peaks  and  valleys  across  the  electrode 
in  the  strongly  coupled  regions.  The  lowest  complex  disper¬ 
sion  curves  for  PZT-7A  plates  have  been  determined  and  will 
be  employed  in  the  analysis  of  energy  trapping  in  the  lower 
frequency  range,  which  has  been  observed  experimentally.  An 
analysis  of  PZT-7A  plates  driven  into  thickness-extensional 
trapped  energy  vibrations  by  the  application  of  a  voltage 
across  strip  electrodes  is  being  performed.  The  plate  waves 
in  the  unelectroded  region  that  carry  energy  away  from  the 
trapped  mode  are  being  included  in  the  treatment.  All  pre¬ 
vious  treatments  were  limited  to  the  consideration  of  an 
eigenvalue  problem;  the  current  effort  will  enable  consider¬ 
ation  of  the  inhomogeneous  forced  vibration  problem. 
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PHYSICAL  PROPERTIES  OF  EPOXY  RESINS  AND  COMPOSITES 

Senior  Investigator:  S.  S.  Sternstein 

This  project  is  concerned  with  the  mechanical  charac¬ 
terization  of  epoxy  neat  resins,  including  particularly  the 
effects  of  moisture,  temperature  and  stress  histories  and 
ultimately  the  dependence  of  the  properties  of  high  perform¬ 
ance  composites  on  neat  resin  behavior . 

The  first  phase  of  this  work  has  involved  the  applica¬ 
tion  of  inhomogeneous  swelling  theory  developed  previously 
for  filled  elastomer  systems  to  the  water  swelling  of  a  neat 
resin  containing  isolated  carbon  fibers.  A  principal  problem 
area  has  been  the  development  of  a  suitable  consitutive  equa¬ 
tion  for  describing  the  thermodynamic  interactions  of  water 
with  neat  resin  and  the  subsequent  mechanical  properties  of 
the  swollen  resin.  Extensive  data  from  the  literature  have 
been  used  to  develop  a  model  constitutive  equation.  Such  an 
equation  has  now  been  formulated,  and  it  appears  to  fit  the 
available  data.  Briefly,  the  model  employs  a  standard  lat¬ 
tice  model  calculation  for  the  combinatorial  entropy  of  mix¬ 
ing  and  a  concentration-dependent  Flory-Huggins  parameter 
for  the  excess  free  energy  of  mixing . 

The  mechanical  work  function  required  for  the  inhomoge¬ 
neous  swelling  field  problem  has  been  developed  from  the 
Gaussian  theory  of  rubber  elasticity  and  utilizes  a  concen¬ 
tration-dependent  crosslink  density.  In  effect,  this  allows 
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the  modulus  to  depend  explicitly  on  water  concentration. 

This  is  to  be  contrasted  with  the  usual  treatment  of  rubber 
swelling  in  which  modulus  and  crosslink  density  do  not  de¬ 
pend  explicitly  on  concentration  but  rather  implicitly 
through  the  volumetric  swelling  ratio  (a  simple  dilution 
calculation) .  It  is  not  suggested  that  the  revised  Gaussian 
model  has  any  particular  molecular  significance;  however,  it 
fits  the  data  well.  One  can  conjecture  that  the  model  is 
that  of  a  two-crosslink  density  network,  one  being  due  to 
the  true  chemical  crosslinks  in  the  epoxy  resin  and  the  other 
a  physical  crosslink  network,  e.g. ,  due  to  hydrogen  bonds. 

It  is  the  latter  network  which  is  destroyed  by  water  absorp¬ 
tion  and  results  in  reduced  modulus. 

The  inhomogeneous  swelling  theory  has  been  revised  for 
application  to  cylindrical  inclusions,  e.g.,  carbon  fibers 
and  sample  calculations  have  been  performed  using  the  new 
constitutive  equations.  It  is  found  that  the  overall  swell¬ 
ing  ratio  of  the  epoxy  is  little  affected  by  the  presence  of 
a  carbon  fiber.  However,  the  local  stresses  and  strains  in 
the  vicinity  of  the  interface  (between  the  fiber  and  epoxy 
matrix)  are  strongly  influenced.  Thus,  a  global  view  of 
swelling  in  composites,  e.g.,  by  weight  gain  studies,  would 
appear  as  a  simple  reduction  of  water  absorption  in  propor¬ 
tion  to  the  fiber  content.  However,  a  local  view  leads  to  a 
much  deeper  insight  into  the  role  which  water  absorption 
plays.  For  example,  it  is  found  that  the  radial  stress  at 
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the  fiber  interface  generated  by  inhomogeneous  swelling  can 
be  as  high  as  50%  of  the  epoxy  tensile  strength. 

The  calculations  have  also  shown  that  under  appropriate 
conditions  determined  by  the  constitutive  parameters  the 
water  concentration  at  the  interface  is  higher  than  at  further 
distance  into  the  matrix.  This  is  in  sharp  contrast  to  swell¬ 
ing  in  a  filled  rubber  where  the  swelling  at  the  interface  is 
markedly  reduced. 

Additional  calculations  on  swelling  instability  (e.g., 
void  formation)  are  in  progress.  A  publication  on  the  swell¬ 
ing  field  calculations  is  being  prepared. 

The  second  phase  of  this  work  involves  the  dynamic  me¬ 
chanical  characterization  of  swollen  neat  resins.  Briefly, 
the  water -equilibrated  samples  are  being  subjected  to  rapid 
temperature  scans  and  the  dynamic  mechanical  loss  factor  is 
being  obtained  at  frequencies  from  1  Hz  to  100  Hz.  These 
data  will  be  required  for  further  analysis  of  the  water  swell¬ 
ing  of  epoxy  resins. 
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FATIGUE  IN  COMPOSITE  t4ATERIALS 

Senior  Investigators:  E.  Krernpl 


During  this  reporting  period  activities  were  mostly  de¬ 
voted  to  the  design  of  a  biaxial  fatigue  test  specimen  and 
associated  fixtures.  The  Cyclic  Strain  Laboratory  is  equipped 
with  an  MTS  tension-torsion  hydraulic  test  system,  and  we 
plan  to  utilize  this  testing  machine  in  our  experiments. 

Much  thought  has  been  given  to  means  by  which  we  can  insure 
that  specimens  will  not  fail  at  the  mounting  grips  when  test¬ 
ed  with  this  equipment. 

The  proposed  specimen  is  shown  in  Figure  31.  It  is  made 
by  wrapping  graphite-epoxy  tape  at  ±45°  around  an  aluminum 
mandrel.  At  present  four  layers  of  tape  are  used.  At  each 
end,  90°  layers  are  built  up  and  subsequently  machined  to 

the  taper  shown  in  Figure  31. 

The  specimen  is  made  to  fit  into  the  fixture  shown  in 
Figure  32.  The  fixture  includes  a  circular  main  part  (1) 
which  bolts  to  the  testing  machine  crosshead.  A  circular 
cover  plate  (2)  wedges  the  specimen  into  the  main  part  through 
the  circular  plug  (3)  with  square  ends  whicn  fits  into  a  cor 
responding  opening  in  the  main  part  (1)  and  by  the  split 
wedge  (5)  and  the  ring  wedge  (4) . 

The  specimen  is  mounted  by  tightening  the  bolts  on  the 
circular  cover  plate  (2) .  Then  hydraulic  pressure  is  ad¬ 
mitted  through  a  separate  part  which  pressurizes  the  plug 
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SPECIMEN 
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snd.  kssps  thG  w©dgs  ssssitibly  undsjr  const3.nt  piTsssuirs.  It 
pjTGvsnts  ths  loosBnin^  of  ths  fit  du,s  to  3rslGX3.tion  in  ths 
specimen. 

The  axial  force  is  transmitted  through  the  wedges. 
Torque  transmission  relies  entirely  on  friction  between  the 
plug  (13)  and  the  specimen  and  the  split  wedge  (5)  and  the 
specimen. 

The  plug  (3) -has  a  central  bore  which  can  be  drilled 
through  the  square  ends.  This  bore  will  permit  the  addition 
of  internal  pressure  as  an  independent  load  system  in  addi¬ 
tion  to  the  axial  loading  and  the  torque  presently  planned. 

Prototype  test  specimens  have  been  made  with  no  diffi¬ 
culty,  and  machining  of  the  fixture  is  nearing  completion. 
Preliminary  tests  are  planned  with  the  MTS  tension-torsion 
system  in  the  next  reporting  period. 

In  parallel  with  the  design  and  fabrication  work  de¬ 
scribed  above,  uniaxial  fatigue  tests  on  { 0/±452/90}2s  and 
{90/+452/0}2s  graphite  epoxy  12-ply  laminates  as  reported 
during  the  last  period  continued,  with  special  consideration 
of  frequency  effects. 
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TRANSVERSE  THERMAL  EXPANSION  OF  CAFIBON  FIBER/EPOXY  COMPOSITES 

Senior  Investigators:  R.  J.  Diefendorf 

C.  LeMaistre 


It  is  important  to  be  able  to  predict  physical  proper¬ 
ties  of  laminates  from  constituent  properties  by  means  of 
micromechanical  analysis.  Several  authors  have  developed 
micromechanical  theories  and  have  been  able  to  obtain  reason¬ 
able  correlation  with  experimental  data.  This  is  particular¬ 
ly  true  of  elastic  and  thermal  properties  with  isotropic  con¬ 
stituents  .  Moreover,  the  longitudinal  TCE,  de¬ 

rived  by  Greszczuk^^^  and  Schapery^”^^  is  even  effective  for 
composite  systems  which  have  anisotropic  fibers  (graphite/ 
epoxy).  Their  expression  for  a,,  is: 


^_mmmrrr  n) 

“11  '  E„(l  -  Vj)  +  EjVj 

where  a,  V,  E  stand  for  TCE,  volume  ratio  and  Young's  modulus; 
and  m  and  f  subscripts  stand  for  matrix  and  fiber,  respective¬ 
ly.  This  equation  assumes  that  Poisson's  ratio  of  the  fiber 
and  matrix  are  equal.  The  good  agreement  is  obtained  because 
the  fiber  properties  dominate  in  the  0°  orientation,  and 
fiber  properties  are  measured  parallel  to  the  fiber  axis. 

Prediction  formulas  for  transverse  TCE  do  not  fare  as 
well  for  graphite  reinforced  composites  as  they  do  for  other 
composites  due  to  the  anisotropy  of  the  fibers.  The  theories 


*Nos.  in  raised  parentheses  refer  to  the  References  listed  on 
pages  98  and  99. 
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to  be  presented  in  the  following  paragraphs  make  the  follow¬ 
ing  basic  assumptions; 

1)  The  ply  is  macroscopically  homogeneous,  linearly 
elastic  and  generally  orthotropic » 

2)  The  fibers  and  matrix  are  linearly  elastic,  homo¬ 
geneous  and  isotropic . 

3)  There  is  complete  bonding  at  the  interface  of  the 
constituents,  and  there  is  no  transition  region 
between  them. 

4)  The  ply  is  initially  in  a  stress— free  state. 

(  g  ) 

5)  The  fibers  are  regularly  spaced  and  aligned. 

1  (6) 

A  mechanics  of  materials  approach  was  used  by  Greszczuk 
to  determine  and  a22-  The  quantity  was  presented  in 

Equation  1  and  a22  is  expressed  as  follows: 


a 


E  B  +  a  Ed  -  B)  ] 

22  £200^  mm 

-  23)  +  2a^B  -  -  aj  (1  -  23) 


(2) 

(3) 


_ _ x:f _ 

^O  [E^d  -  23)  +  2E^B] 

E„  =  E^B  +  E^d  -  3) 

2  0  m 


This  formula  assumes  circular  filaments.  Schapery  devel¬ 
oped  bounds  on  the  TCE  of  isotropic  constituent  composites  by 
employing  extremum  principles  of  thermoelasticity.  As  shown 
before,  Equation  1,  Schapery 's  formula  for  is  the  same  as 

that  derived  by  Greszczuk.  Schapery  also  presents  a  simple 
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expression  as  an  upper  bound  for  ^22  Wang  Fo- 

®22  =  *  Wm  *  Vm' 

The  results  from  formula  (7)  compare  well  with  experimental 
results  from  fiber  glass  composites. 

Fahmy  and  Ragia-Ellozy and  Kahn^^^^  applied  finite 
element  analysis  to  determine  and  ^^2  ^  lamina.  This 

approach  provides  a  rigorous  solution  for  ^22'  requires  a 

sophisticated  microstructural  analysis  of  the  composite  sys¬ 
tem  and  does  not  provide  a  readily  accessible  equation  for 
determining  0.22  fj^om  constituent  properties.  Kahn  also  de¬ 
velops  an  equation  for  0.22'  treating  the  fiber  as  a  solid  rod 
of  radius  'a'  and  the  matrix  as  a  tube  of  inner  radius  'a' 
and  outer  radius  'b'.  The  matrix  tube  is  subjected  to  an  in¬ 
ternal  pressure,  P,  whose  sign  depends  on  the  thermal  mis¬ 
match  between  fiber  and  matrix; 


^22  “m  AT 


(£  )  =  A.[  2PJL.  _  y  ^  ] 

^^em^b  E  ^1  -  X  zm 

m 


V, 


E^d  -  V^) 

+  Ed-  V.)  (1  +  V.)  +  E  V  (1  +  V  )^1  -  V 
m  mf  rrrm  r 


(8) 

(9) 

(10) 


These  equations  result  in  a  curve  which  falls  between  the  rule 
of  mixture  lower  bound  solution  and  Schapery's  upper  bound 
equation.  Finally,  Levin ^  provides  an  anisotropic  analysis 
for  the  thermal  expansion  coefficients.  No  comparison  is  made 
to  experimental  results  for  almost  all  analyses. 
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In  an  attempt  to  evaluate  the  applicability  of  these 
theories,  thermal  expansion  specimens  were  fabricated  and 
tested  on  a  push-rod  dilatometerc  The  material  selected 
for  this  investigation  was  the  Fiberite  1048  AE  graphite/ 
epoxy  system.  Components  consist  of  Union  Carbide  T-300 
graphite  fiber  and  Fiberite 's  121°C  (250°F)  48A  resin  in 
prepreg  form.  Test  specimens  were  fabricated  on  a  Dake  75- 
ton  press  with  19"  x  19"  programmable  heated  platens.  The 
cure  cycles  were  programmed  with  a  Data  Trac  and  fully  moni¬ 
tored.  Nine  specimens,  6"  x  6",  were  fabricated  with  vary¬ 
ing  fiber  volume,  from  which  five  expansion  samples,  2"  long, 
approximately  1/2"  in  highth  and  thickness  varying  depending 
on  fiber  volume,  were  cut.  All  nine  specimens  contained  32 
plies  unidirectional  prepreg  cured  at  various  pressures, 
quantity  of  bleeders  and  excess  resin  to  achieve  the  required 
resin  content.  Neat  resin  was  cured  in  the  same  manner  as 
the  laminates  under  vacuum  pressure  to  remove  porosity.  Sim¬ 
ilar  dimensions  were  used  for  the  resin  TCE  specimens.  Dog- 
bone  tensile  specimens  were  fabricated  for  determination  of 
tensile  modulus  and  Poisson's  ratio  of  the  resin. 

TCE  specimens  were  conditioned  in  an  air  circulating 
oven  at  110®C  for  24  hours  and  stored  in  a  desiccator  until 
tested  in  an  attempt  to  eliminate  moisture  induced  variations 
in  the  thermal  expansion  measurement Specimens  were 
v;eighed  before  and  after  testing  to  determine  moisture  loss 
during  testing. 
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Thermal  expansion  measurements  were  performed  on  a  modi¬ 
fied  Thermo-Physics  Corporation  push  rod  dilatometer.  Sam¬ 
ples  were  heated  from  room  temperature  to  lOO^C  at  a  low  rate 
of  0.75®C/min.  to  minimize  discrepancies  which  may  arise  from 
the  variation  in  thickness  and  thermal  conductivity  of  the 
different  fiber  volume  samples.  The  linear  TCE  was  derived 
from  a  continuous  recording  of  length,  utilizing  a  precision 
LVDT,  with  resulting  accuracy  of  0.1  x  10  /°C. 

Two  methods  were  used  to  determine  fiber  volume?  1)  vis¬ 
ual  count  method  using  exploded  photo-micrographs  of  the  com¬ 
posite  microstructure  and  2)  Boeing's  MS8-212A  nitric  acid 
digestion  method.  Densities  were  determined  utilizing  ASTM 
0792  "Tests  for  specific  gravity  and  density  of  plastics  by 
displacement.  Method  A". 

Mechanical  testing  was  performed  on  an  Instron  screw- 

type  universal  testing  machine  at  a  crosshead  rate  of  0.005"/ 

min.,  and  strain  data  was  recorded  on  a  Data  Logger.  The 

1048  AE  composite's  mechanical  property  data  used  were  re- 

(14) 

suits  of  early  testing  by  Kowalski 

The  results  of  individual  testing  indicate  that  the  24 
hour  thermal  conditioning  of  the  TCE  specimens  was  effective 
in  minimizing  the  moisture  response  on  expansion.  The  typi¬ 
cal  peaking,  probably  caused  by  the  T  depression,  evident 
in  past  work  (^2,13,15)  absorbed  moisture,  was  not 

noted.  Weight  loss  related  to  moisture  desorption  during 
testing  was  less  than  .03%  in, all  samples  tested. 


92 


Table  XII  shows  average  transverse  -TCE  from  30°C  to  100 °C 
versus  fiber  volume.  These  results  are  plotted  on  Figure  33. 
From  these  data  points,  an  estimate  of  ^22  individual 

fiber  was  calculated  by  utilizing  a  least-mean-squares  extrap¬ 
olation  to  =  100%.  It  is  noted  that  a  modified  rule  of 
mixtures  behavior  is  assumed  and  that  a  straight  line  is  valid 

to  approximate  ^22'  this  analysis  it  is  estimated  that 

—  6 

6^22  T-300  graphite  fiber  has  a  value  of  8.8  x  10  /®C  in  a 

temperature  range  of  30°C  -  100°C.  This  result  is  similar  to 

a  result  obtained  by  Ishikawa^  '  who  estimates  a  value  of  8 . 6 

X  10”^/®C.  However,  it  is  important  to  note  that  Ishikav/a's 

method  appears  to  be  self-fulfilling  in  that  0-22  determined 

by  extrapolating  values  from  a  theoretical  curve  and  plugging 

the  results  back  into  the  theory  to  prove  its  fit. 

In  order  to  correlate  this  estimated  value  of  0.22 

300  to  Schapery's  upper  bound  solution,  elastic  constants  of 

the  matrix  are  required.  Test  results  indicate  that  = 

55  X  10 ^  psi  and  v  =  .35.  These  values  agree  with  typical 
^  in 

^  (16) 

values  of  E  and  for  other  epoxy  resins  used  in  composites 
m 

The  upper  bound  equation  presented  by  Schapery  was  modi¬ 
fied 


a 


22 


(1  - 


V  )  a  V  + 
m  mm 


(1  +  V 


f  )  '^f  '^f  ~ 

^12^  22  ^ 


a,  (V.  V-  +  V  V  ) 
3.  ^  mm 


to  include  fiber  anisotropy,  a,  is  defined  by  Equation  (1) 

and  would  appear  to  be  close  to  the  value  high  modulus 

^11 

fibers  at  moderate  or  high  volume  fractions.  Direct  applica¬ 
tion  of  this  equation  gives . a  theoretical  curve  that  is  higher 
than  all  experimental  values,  because  the  boundary  condition 
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TABLE  XII 


Specimen 

Number 

TEST  RESULTS: 

TRANSVERSE 

TCE  VERSUS  FIBER 

VOLUME 

Density 

(g/cc) 

(nitric  digest) 

(micro-photo) 

TCE 

(ave. 

(xlO  V°C) 
:30-100°C) 

11 

1.542 

67.2 

23.5 

12 

1.595 

70.1 

27.2 

13 

1.454 

— 

25.1 

14 

1.586 

64.2 

26.6 

15 

1.585 

69.2 

25.2 

ave. 

67.7 

68.7 

25.5 

21 

1.553 

30.4 

22 

1.563 

65.1 

29.7 

23 

1.524 

64.7 

29.7 

24 

1.533 

64.6 

28.2 

25 

1.580 

69.3 

26.8 

ave. 

65.9 

63.0 

28.9 

31 

1.628 

67.3 

31.9 

32 

1.550 

63.9 

28.7 

33 

1.543 

66.1 

30.1 

34 

1.514 

62.4 

33.7 

35 

1.542 

63.8 

31.0 

ave . 

64.7 

64.2 

31.1 

41 

1.530 

61.9 

32.3 

42 

1.557 

62.1 

32.9 

43 

1.579 

62.9 

29.7 

44 

1.557 

62.4 

33.9 

45 

1.530 

63.8 

31.3 

ave. 

62.3 

60.9 

32.0 

51 

1.532 

59.2 

36.2 

52 

1.537 

59.4 

31.7 

53 

1.525 

58.8 

35.4 

54 

1.500 

57.4 

34.7 

55 

1.513 

58.1 

33.4 

ave. 


58.6 


56.3 


34.3 


94 


TABIiE  XII  continued 


Specimen  Density  f 

Number  (g/cc)  (nitric  digest) 


61 

1.470 

4^ 

o 

62 

1.451 

50.9 

63 

1.443 

50.4 

64 

1.454 

51.2 

65 

1.449 

51.3 

ave. 

49.4 

71 

1.501 

52.1 

72 

1.492 

54.4 

73 

1.489 

54.4 

74 

1.489 

54.1 

75 

1.478 

54.9 

ave. 

54.9 

81 

1.354 

46.9 

82 

1.401 

46.8 

83 

1.407 

— 

84 

1.417 

1 — 1 

00 

85 

1.431 

52.0 

ave. 

48.4 

91 

1.393 

40.4 

92 

1.343 

38.4 

93 

1.328 

37.3 

94 

1.342 

38.7 

95 

1.389 

41.3 

ave. 

39.2 

101 

1.185 

102 

1.201 

103 

1.200 

00.0 

104 

1.213 

105 

1.205 

TCE  (xlO^V^C) 
(micro-photo)  (ave  o  g  3Q-1Q0^C) 

34.7 

37.1 
39.9 
40.5 

39.1 

51.5  38.2 


38.3 

36.6 

38.1 

39.1 
35.0 

53.6  37.4 


40.0 

42.4 

38.6 

41.2 

38.8 

49.6  40.2 


46.1 

44.1 
41.0 
44.4 
43.3 

36.6  43.8 


51.9 

53.9 
52.5 

52.8 

54.8 


ave. 


53.2 


CE  versus  Fiber  Volume 
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at  Vj  =  1  is  not  satisfied.  In  order  to  satisfy  this  condi¬ 
tion,  the  following  must  be  trues 


lim  = 

V^-^1 

(1  -I-  V-  )  a-  -  a 

£12  £22  ^  ^12 

(11) 

1,  ^22 

(X  £  —  Ct^  t  ^-F  ” 

^22  22  12  22  12 

(12) 

®a  "  “f.. 

(13) 

The  resulting  equation  is: 

a.-  =  (1  -  v^)a  V  +  (V.  -  V  V  )a. 
22  mmm  £  mmr 


(7b) 


22 


Input  values  used  in  the  theoretical  prediction  plotted  on 
Figure  33  are: 


.35 

.30 

(Ref.  15) 

V  = 

m 

^12 

m 

53.2  X  10' 

■^/°c 

Uf  = 

^22 

8.8 

X  10“^/°C 

E  = 
m 

.55  X  10^ 

psi 

^f  = 

32. 

5  X  10^  psi 

It  is  seen  from  Figure  33  that  excellent  correlation  exists 
between  experimental  data  and  the  results  calculated  from 
Equation  (7b) .  The  reason  for  the  need  to  modify  the  form 
of  the  equation  falls  back  on  the  hidden  assumption  in  the 
derivation  of  Equation  (7a)  in  the  third  term  that  the  con¬ 
stituents  are  isotropic  in  nature.  Equation  (7a)  assumes 

K 

volume  ratio  of  the  bulk  modulus  (i.e. ,  liK  V  )  and  K  is  de- 

n  n 


that  a  =  St T  from  isotropy  and  a  -  where  K  is  the 

a  IX  ^  K 

de 

fined  to  be; - - - — 


^11  ^  ^22  ^33  "  “v^ 


Ka  is  defined  in  the  same  manner.  In  Equation  (14) ,  is 
equal  to  +  ^22  ^33'  which  in  the  case  of  isotropic 
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constituents  can  be  reduced  to  Substituting  this 

isotropic  constituent  condition  in  Equation  (14)  and  assuming 

that  ^  ~  Equation  (7)  is  derived  by  energy  principles. 

■  E  K 
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